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16000 PILE CYCLIC LOADS AND LOAD RATES 
| KEY WORDS: Capacity; Clays; Compression; Cyclic loads; Loading rate; yay hs 
| Pile load tests (cyclic loading); Soils; Steel piles; Stratigraphy; Tension; = 
odi bas soneduseib Yo od) noisibbs ai hoa beteslave 
| ABSTRACT: Two series of axial load tests were performed on four 14-in. (356-mm) 
| diam, open-end, steel pipe piles at an interval of about 320 days. Pile lengths of 40 ft a 
or 50 ft (12.2 or 15.2 m) were installed below conductors driven to depths ranging | Ny 
_ from 115 ft to 320 ft (35 m to 98 m) into a strong underconsolidated clay. Each pile 
- was subjected to as many as 26 tests. Data presented includes: (1) Compression and 
_ tension tests; (2) tests performed at different times after driving and after previous 
{| tests; and (3) incremental loading and constant rate of loading. Pile capacity increased 
+ 40% to 75% when the loading rate increased by three orders of magnitude. The one- 
- way cyclic loading applied in this study did not effect the ultimate capacity, but large 
displacement began to — when the maximum cyclic load reached 80% of the 
| REFERENCE: Kraft, Leland M., Cox, William R., and Verner, Edward A., 
Tests: Loads and ‘Load Rates,” 


15980 LONGITUDINAL VIBRATIONS OF EMBANKMENT DAMS > 


+5 KEY woans. Earth Earthquakes; Embankments; Soil dynamics 


y ABSTRACT: This paper presents an analytical method to study free and forced 
y > longitudinal vibrations of embankment dams. Both shear and dilatational deformations 
: | are taken into account, and the dam is modeled as a linear homogeneous prism with a 
= _ wedge-shape cross section, bounded by two vertical (abutment) and one horizontal _ 
Pe (riverbed) planes. Numerical results demonstrate the effect of the length and height of 
' _ the dam on its. natural frequencies and its modal displacement and strain shapes. It is _ 
hi shown that shear deformations are more important in relatively long dams, whereas the 
' opposite is true with dams built in narrow canyons. The method is evaluted through | 
ig | _ two case histories involving an earthfill and a rockfill dam. Close agreement exists 


accelerations within the dam is badly predicted. 

At nominal :neresate in cost by testes cifleren 

_ REFERENCE: Gazetas, George, “Longitudinal Vibrations of Embankment Dams,’ 

_ Journal of the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT1, Proc. 

of the performance of piles subjected axial Dante are fer 
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15979 CYCLIC AXIAL ‘RESPONSE OF SINGLE FILE 

KEY WORDS: Bearing capacity; Consolidation; Cyclic loads; Foundations; 
Offshore structures; Pile settlement; Piles Shear strain; Skin 


friction; Soil mechanics 
_ ABSTRACT: A number of solutions have been obtained for a typical offshore pile, to 
determine the influence of a number of the input parameters on the computered cyclic 
response. More significant parameters include the critical shear strain at which < 
significant degradation of skin friction occurs, the rate of loading on the pile, and the 
‘ distributions of static skin friction and soil modulus along the pile. The analysis 
predicts a gradual decrease in cyclic stiffness of the pile with increasing numbers of 
cycles and increasing cyclic load level, but a very sudden decrease in ultimate load 
capacity once the (half-peak-to-peak) cyclic load level exceeds 50% to 60% of the 
1 ultimate static load capacity. These theoretical findings are broadly confirmed by the 
REFERENCE: Poulos, Harry G., “Cyclic Axial Response of Single File,” Journal of 
_ the Geotechnical Engineering Division, ASCE, Vol. 107, No. GTI, Proc. Paper 
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15985 ANCHOR RESPONSE TO UPLIFT LOADING 

KEY WORDS: Anchors; Foundations; Ocean Soil mechanics 
ABSTRACT: A number of factors which influence soil- anchor interaction are 
evaluated, and in addition, the extent of soil disturbance and the influence of pore 1am : 
water pressures during testing are assessed. Response under repeated and sustained- | 
“repeated loading is shown to be represented by a family of hyperbolic curves, and a 
failure criterion is determined by the critical | strain condition below which a 
satisfactory response can still be subsequent to repeated loading. A strain dependent 7 
mechanism is used for defining behaviour under varying loading conditions which 

- permits an assessment of anchor life when subjected to a package of projectd — 
environmental conditions. Where possible the results are compared to those of other — 

workers in the anchor field, and ee — — may lead to a Setter es 


» 


: REFERENCE: Andreadis, Andress, Harvey, Roger C., and Burley, Eldon, ae 
Anchor to uplift Loading,” Journal of the Geotechnical Engineering 
pp. 59-78 


15992 CONSTRUCTION OF CANAL ON COLLAPSING SOILS 0.2 nee 
KEY WORDS: Canals; Collapsible soils; Depth; Field density; Liquid limit; 
Pneumatic tired rollers; Pondweeds; Soil consolidation tests; Subsidence __ 
| apsTRACT: Soil subsidence may occur because of withdrawal of fluids (usually deep | 
t subsidence), the drainage of peat lands, or by the application of water to moisure- 
ok, deficient, low-density soils (usually near-surface subsidence). Only near-surface 
' subsidence is analyzed in connection with the investigations for, and construction of, _ 
_ the San Luis Canal where parts of the alinement passed through areas of low-density, — 
-moisture-deficient soils. Ponding of the critical subsidence areas delineated by studies 
. — an effective methods of collapsing and densifying the soils prior to construction of 
sca canal. The canal invert and the embankment were densified with a heavy : 
_ pneumatic roller during the construction phase. The San Luis Canal was constructed Xe 
during 1963 to 1968 and has functioned with minor to date (1979). 


| REFEREN 
Journal of the Geotechnical Engineering ASCE, Vol, 107, No. GT1, Proc. 
| Paper 15992, January, 1981, pp. 79-94 


15998 COMPUTER SIMULATION OF ‘CREEP OF PCLAY 

KEY WORDS: Clays; Computerized simulation; Creep; Creep rate; Energy; ‘v 
4 bands; Illite; Microstructure; Soil mechanics; Soil wae Stochastic > 


ABSTRACT: The variable character o of the: particle” contacts in illitic 
- implies a spectrum of bond strengths. When a sufficiently high external shear stress a 
applied some weaker interlinking particle groups are transformed to “domains” which 
act as slip units. In creep induced by a constant deviatoric stress, the number and - 
heights of thermally activable barriers to slip will vary in space and time. This is 

_ encompassed by a stochastic model, a simple version of which was used in a computer 

_ simulation. The dependence of the strain on time was found to be logarithmic, in os 

_ REFERENCE: Pusch, Roland, and Feltham, Paul, “Computer Simulation of cal of | 

_ Clay,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 107, No. GTI, 
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cman with h 13 other participants to provide data on the axial-load capacity 
of deep foundation piling in strong clays that are underconsolidated. Two series 
_of axial load tests were performed in 1975 on four 14-in. (356-mm)-diam, open-end, _ 

steel pipe piles with an interval of some 320 days between each series. The | 7 


7 site was located one mile south of Empire, La. (Fig. 1), ie 


‘The scope of the initial testing program was limited to pile capacities eter 
_ load conditions. Additional information about pile performance could be 
_ obtained at a nominal increase in cost by including tests at different loading 
rates, so the project scope was poreareagieagen i the project. The data from 
the two test series offer an opportunity to gain an improved enderstending 
of the performance of piles subjected to axial loads. Data are available for: 
(1) Compression and tension tests; (2) tests performed at different times after 
- driving and after previous tests; (3) incremental loading and constant | rate of — 4 
— loading tests; (4) tests loaded to failure in periods ranging from a few seconds 
to hours; and (5) cyclic load tests with periods typical of ocean wave loading. | 
_ Pile performance was studied in terms of: (1) Load-deformation (pile-soil stiffness) 
_ response at the top of the test pile; (2) cumulative pile deformation; (3) —— 
of load transfer along the pile with load application; (4) redistribution of load 


a. Megr., Special Projects Group, McClelland Engrs., , Inc., 6100 Hillcroft, Houston, Tex. se 
“Vice and Mgr., Offshore Engineering, McClelland Engrs., Inc., 6100 Hillcroft, 
Pres. and Mgr., Offshore Technology, TERA, Inc., Houston, Tex. 
~ Note.—Discussion open until June 1, 1981. To extend the closing date one month, 
-awritten request must be filed with the Manager of Technical and Professional Publications, _ 
= Manuscript was submitted for review for possible publication on March 25, 1980. : 
- This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
| the American Society of Civil Engineers, © ASCE, Vol. 107, Aptos GTI  Janeery, 
ISSN 0093-6405 /81 /0001-0001 /$01.00. 
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FIG. ¥ 1.—Test Site and Boring Plan (1 ft = 0.305 m; iin mile = 
transfer along the pile with cyclic loading and unloading; and (5) ultimate ‘elle 
‘This paper describes the soil conditions, the piles, the test program, the reduced _ 
load test data, and an evaluation of these 


; Soil Borings.—The soils investigation included: (1) Three undisturbed sample | 
q “borings a at b the § site to determine soil ‘Stratigraphy and to ) obtain samples for 


a a boring | to additional definition of soil (3) field and laboratory ¥ 


= on recovered samples to determine pertinent physical properties of the — 
soils; and (4) engineering analyses of the field and laboratory data for fener 


The soils investigation also included in situ remote (6,7) tests. 


z= of the test program and to develop pile capacity predictions. 


of classification tests and strength 1 measurements from 
| unconfined compression tests, unconsolidated-undrained tests, and miniature © 
= tests were performed. In addition, stress-strain curves and measured shear | iA 
trengths were determined from six triaxial tests, static ‘simple 


| 
al 


PILE LOAD TESTS 


shear test tests, five cyclic | simple : shear tests, three — tests, 16 friction 
tests between clay and steel (eight of these tests were on remolded specimens). ¥ 


7 Pa _ Stratigraphy.—The borings drilled at the test site revealed a stratigraphic | 


equence containing two major clay strata in which the pile tests were performed. - 


_ These strata are a firm, highly plastic, uniform clay extending from 76-ft to 
178- ft (23-m m to. 54-m) penetration; and a stiff to very stiff, ; gray, highly a 


the shear strength profile are summarized in 2. Four 40-ft 
to 50-ft thick zones of clay between 115-ft and about 360-ft penetration were 7 

from data it in which t to the pile tests. These test zones 

sah 


MEAN REGRESSION 
STRENGTH 
a INTERPRETED MEAN 
FROM ZON 
O INTERPRETED MEAN STRENGTH 


INTERPRETED 
STRENGTH 
PROFILE 


REGRESSION 

LINETO 


sainshditintts at Empire Te: _ st Site at ft = = 0: 305 m; 1k ksf = 47.91 kN /m? 4 


are denoted by “ “Pile 1,” “Pile 2,” **Pile 3, and “Pile 4 i in Fig. Table. 


4 
—_ ‘1 summarizes the range in interpreted soil properties for each test zone. mor’ 


‘oman in Zone | is 1.20 ksf. Bas on) al 
The material in Zone 2 is more variable, based on resistivity, gat gamma ray, 
and sample logs, than the soil in Zone 1. Many thin seams of sand and silt 


were revealed in the clay. The mean shear strength from the interpreted profile 


resistivity, and gamma- ray logs. Sample recovery was “poor, , due 


= to its layered condition, and strength data are concentrated at the | 


> bottom portion of this zone. Test results from 205-ft, 302-ft, 305-ft, 310-ft, 
and 320-ft penetrations were larger results at other ‘penetrations. These 


. a Resistivity and gamma-ray logs, : sample descriptions, and Atterberg | limit data | 


JANUARY 198 


larger re results indicate a mean strength. of about 2.6 ksf. Judgment interpretation, 
_ relying on miniature vane data largely in a zone at 250-ft penetration and another ee 


from to 345- ft penetration, suggests a mean | strength of about 2. 


The resistivity ‘and gamma- ray logs indicate sandy material 346-ft to 
352. ft penetration in Zone 4. Sample descriptions show that the clay in this 4 
zonei is layered with sand and silt seams. The strength data exhibited considerable 
scatter in Zone 4 due to variable sample disturbance caused by layering and Jy 
in physical Properties of the material. A weighted mean strength 


- ierpretation , based largely ¢ on n data at 250 ft and from 340 ft to 345 ft, — 


1 of Soil Conditicns at 
Zoe | | Zoe 
Parameters” 
_ 


_ Mean undrained shear strength, s,,,, in rh 
a kips per square foot = 
Mean effective overburden pressure, 
in kips per square foot 
angle, ®, in degrees ad 
Earth ‘at rest,” 


| 


- Rebound to compression index ratio, 0.2 
4 


Note: ksf=47.9kKN/m 


a mean value for this test zone as much as 3.95 ksf. An interpreted strength 
_of 3.58 ksf for Zone 4 was used from the profile shown in Fig. 2. euvh 


Test Pues ano Test Sequence 


wen driven to the top of the test zones and then the soil oy in the conductors — 
were drilled out. The four 14-in. (356-mm)-diam, open-end piles were then 
driven into the test zones below the conducters. made on 
ii first series of compute and tension ania load tests were made on 


consisted generally of slow static (incremental loading) tests in compress 

and tension modes. A. few ‘‘quick static’’ tests (constant loading rate) were 
performed selectively on two piles. 
&s The first tests in me second series were made some 320 days after the last 


= 
ma 
6 
0.55 
— 
0.28 | 032 

03 

| 
f —_—_loading for each test. At the end of the static tests in the second series, cyclic 


PILE LOAD T 5 
tests with about a a 15-s period were run in compression mode on each of the - 
_ four piles. At the end of cyclic testing, static tests at varying rates of loading — 
were made. A series of very fast compression loadings were then made under 


_ separate funding provided by Shell Oil ‘Company; these tests will be referred 


shown in Fig. ;& Test information, including the setup time, are tabulated in 
Table 2, which also summarizes the loading rates foreach testtype. 
_ A schematic diagram of the test setup is presented in Fig. 4. The acquired 7 
data consisted of the axial strains at four points along the length of the test : 
section, the pressure in the jacking system, the vertical movement of a weight — 
attached by wire line to a weight at the top of the test section (diaphragm), 7 b 
and the vertical movements of two plates on either side of the top of the 
_‘14-in. diam pipe. These measurements combine to give the loads at four points a 
along the test section (with a backup for the top of the test section), and the _ 


displacement att the top oft the test section with a backup measurement at 


DRIVING DRIVING COMPRESSION | | DRIVING _ 
COMPRESSION | | COMPRESSION TENSION _|_ | COMPRESSION 


TENSION __|| TENSION __ TENSION COMPRESSION 


COMPRESSION | | COMPRESSION COMPRESSION TENSION —_ 


SECOND | | SECOND 
TENSION TENSION _| |COMPRESSION 


CYCLIC 
VARIABLE FAST 
SPEED STATIC] | CYCLIC SPEED STATIC] | CYCLIC | 
1 
[Tests 


FIG. 3.- Sequence 
evel. Py lead wires for the downhole bridges on Pile 2 were broken during 
ah recovery attempts, rendering all bridges for that pile inoperative. =i 
_ The results were sampled with a data logger on a sequential basis during 
testing with an elapsed time between samples of about 1.5 s. The data were 
automatically r soceeeee, and the the tim ne s shift was accounted for in the data reduction 
All displacement measurements were made with transducers srs of a linear variable 
differential transformer type that were individually calibrated. During the first a 
series of static tests, the displacement at the top of the test section was measured ia 
with the backup system, because the primary wire-line system failed to function. — 
These backup measurements were arrived at by measuring the deglactment 
4 at the top of the ‘14 in. . pile and adjusting this value by the calculated amount 


of shortening, or extension, for the long length of pipe between the pile top 

and the diaphragm at the top of the test section. For the longer piles this 

adjustment was of the same order as the top displacement, resulting in considerable 


any errors | in the component measurements - Probable errors 
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PILELOADTESTS 
ons to such amplifications are conservatively estimated to for Pile 
4 +40% for Pile 2, +75% for Pile 3, and +50% for Pile 4. Because of wordy : 


series of tests to result in an error of less than +2%. bus 


x Loads at the top | of the test. ‘section were - measured with a : a strain 4 gage bridge 
at that level and with a backup reading of the jack pressure. Overall accuracy 
_ of the measurement of the load at the top of the test section was estimated, 
| ; == all contributing effects. The probable error in the measured ery 

= the larger of +5 kips or +2 


FIG. of Test Setup 


dur dures, and data ta accuracy were presented | earlier (5) a and | are not t repeated here. : 


Measured AND PREDICTED Put Benavion were unespected. 


‘The data from the two tests series were used to evaluate the effects of rate 
of loading, elapsed time after driving or previous test, and cyclic — 
_ pile capacity and load-displacement response of the pile. Because each pile _ 
was by choice loaded to failure a number of times, the evaluation of the effect — 
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FIRST 
SERIES 


FIRST SECOND 
SERIES SERIES 


1. FOR TEST NO. SEE TABLE 2 
2. MODE OF LOADING § C = COMPRESSION 
3. DASHED LINES REPRESENT YIELD 


Fi FIG. 5. —Chronological Summary of Pile 


TIFFNESS, KIPS / IN. 


st 


FOR TEST NO. SEE TABLE? 
2. MODE OF LOADING : C = COMPRESSION . T = TENSION 
3. “REPRESENT CYCLIC TEST 


Py. 


FIG. 6.—Chronological of Meneured Pile- Soil Stittness 0. 175 
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of parameter subsequent t to the very first t test may biased 

_by the loading history. Therefore, in many cases qualitative rather than quantita- — 
“tive conclusions have to be drawn from the test data. 
ae. _ Comparisons were made between measured and predicted results; the latter 7 

a ; were determined prior to evaluation of the pile test data. No effort was made > 

a in this study to modify any input data to fit computed with measured values 
= to develop a modified predictive method with the advantage of hindsight — q 


‘Historical Summary of Capacity and Pile-Soil Stiffness.—The capac capacity of a 
pile can be defined in terms of either a yield or ultimate load. For this study, 
ultimate load was defined as the maximum load in a test, while yield was defined -: 
na as the load corresponding to the greatest curvature in the load-displacement — 
curve. The ultimate capacities and yield loads of the four piles are presented eS 
in Fig. 5 in chronological sequence of the test events. 


- Pile-soil stiffness values, as observed from the load- -displacement records « of 7 
the tests, vary with the load level at which they are measured and with the | 7 
G definition of stiffness. The pile-soil stiffnesses of the four piles in terms of 
the secant at yield are | in n Fig. 


_ Measured Static Capacity. the peak « of the first 

~ test of each series (Tests _ 1 and 7 in Fig. 5) show that, after adjusting for 
4 end bearing, Piles 2, 3, and 4 of the second series have lower capacities than 
: the first series despite about 320 days of setup for the second series as compared _ 


“to about 7 days for the first series. Only Pile 1 rr a higher capacity during | 
the second series compared to the first series. 
? _ For Pile 4, the first test of both series was a ‘compression test. Although 
a tests should allow a one-to-one comparison of compressive capacity due _ 
= to a large setup time, the amount of end bearing that could be developed in — 
the second series is questionable as the location of the pile tip in the first 
= of | the second series was above its deepest penetration in the first series. 
- 7 summarizes the residual displacements of the pile after each test sequence. ’ 
The capacity of Pile 4 was only 380 kips (1,690 kN) in the first test of the 
f second series at 327 days after the first series compared to 539 kips (2, 400 
g for the very first test. The computed end bearing of only 38 kips (170 _ 
KN) cannot account for the 159 kip (710 kN) lesser capacity in the test with 
_ Acomparison of the first test of the first and second series shows the following 
differences relative to the first series—Pile 1: +21.6%; Pile 2: —10.2%; Pile 
3: —26.8%; Pile 4: —24.1%. To make this comparison with the tension tests Ke 
of the second series a computed end component was deducted from 
ith 
depth and a decrease i in pile length. ry mi dod 
_ The lower capacities during the second series were re unexpected. The cause 
_ for the lesser capacities in the second test series as well as the low capacity 


| 
| 
q secon ests in € Tirst series is probably to residual strengin ¢ 
; Bergemann (2) observed a similar reduction in capacity after an initial load j 
to failure. In this case, the friction capacity was reduced by 33%, and __ 
on was not recovered even after 2 months of setup. beol 


residual” 4 ‘surface, particularly, for Piles 2, 3, 4, was 
created in the ‘first tests. If such a surface was created by the relatively large 
ile movements in the first tests and the reversal in pile movement due to 
ol subsequent tension and compression testing, the second series of tests would - 
have then indicated the “‘frictional’’ resistance developed by 320 ) days of setup 
on this ‘‘prefailed’’ surface. mao fi cae adi ot 
2 _ Some suggestion of residual strength development and possible : strain sensitivity 
was evidenced by the difference in the shape of the post- peak portion of the : 


load applications. Load-displacement curves for all ‘Static tests on each of the 


four piles are presented in Fig. 8. In the first test on each pile, the ex ° 
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the test section exhibited ‘‘strain softening’’ a 
‘These were all compression tests, the degree of softening w was much less | for : 
Pile 1. Subsequent tests in the first series did not reveal any significant strain _ 
Softening, and the load-deflection curves for the first tension and compression | - 
tests in the second series exhibited only a slight strain-softening tendency. — oa 
‘Strain softening behavior in the laboratory usually occurs with the development 3 
of a pronounced failure surface in the specimen, and the shear resistance at 
large strain is a measure of the residual strength. Stress-strain curves from — 
_ simple shear tests show that the susceptibility of the soil to exhibit a strain 
_ softening behavior increased with an increase in sample depth. Similar behavior 
was observed for the direct shear soil-steel adhesion tests. The shape of the — 
Toad- -displacement curves of the pile and the data both 


q 


PILE LOAD — 11 
the contention a “‘residual”’ strength surface developed with the large pile 
_ movements s during the first test. Subsequent tests results, therefore, are a measure — 
5 g of pile response affected by aresidual strength, 
Predicted Static Capacity.—Procedures to compute pile capacities have histori- bi 
ally been developed from data on pile tests loaded to failure slowly (an elapsed 7_ ; 
time of several hours). The first test of both the first and second series was 
_a slow test and, because of the setup time, should be least affected by the 
a loading history. The ultimate axial capacities were | computed with three methods 
to predict shaft friction [American Petroleum Institute (API), Burland, and 
Lambda (A) methods]. For the soil conditions at the site the unit shaft friction, 


4 bpd the API method (1) is equal to the interpreted undrained shear a 


. Burland @ equates the unit } skin friction to an earth pressure | coefficient 
sin in which d= 
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ce. 8.—Load- ‘Displacement Curves (1 = 25.4mm; 1 ‘Vkip = 4.45 KN) 


an is a quasi-effective stress approach in which the average unit | shaft friction 
is related to the undrained shear strength, the effective overburden pressure, 
and an empirical coefficient that varies with pile penetration (10), 
The undrained shear strengths used in the analysis are shown on the interpreted a 
_ profile in Fig. 2. The average angle of internal friction from the laboratory 7 
tests was 25°. Effective overburden pressures equal to the mean of the range 
‘shown i in Table 1 were used for the predictions. These mean effective stresses BI 
,°§ correspond to 90% consolidation for Piles 1 and 2 and 95% consolidation for 
‘The computed results for each of the methods are compared in Fig. 9 with 
- the measured axial capacities from the first compression test in the first series. _ 
‘Comparisons « are not shown for subsequent tests as the measured results were 
| Comparis by the previous loadings to failure that are not accounted for by 7 
- prediction methods. The capacities shown for the \ method are based o 
an for piles i in normally consolidated soils (8). or BENT 
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of the measured load was estimated tuo be kips (23 
C KN). The uncertainty in the computed values of shaft friction for a given method aa 
g can be related to the uncertainty in the undrained shear strength for the API = 


_method, the undrained shear strength and effective stress for the \ method, - 
 ~. the effective stress and angle of internal friction, ®, for Burland’s method. © 
| The uncertainty i in computed end bearing is related to uncertainty in the bearing 
= factor, N., and the undrained shear strength. The value of N. ranges . 
_ variations in computed values of pile capacity for the uncertainties in 
the soil parameters that are reflected by the range in properties in Table 1 
and for independence between uncertainties are about + 10% for the API method - 
for th the piles in the upper two zones and 20% for the other two piles, =o 


MEASURED CAPACIT 
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PREDICTED CAPACITY, Kips 
: FIG. 9. aniaaapeidaiaeh of Measured and Predicted Ultimate Capacities for First Com-_ 
pression Tests of First Test Series (1 kip = 4.45 kN) 
+10% the method, and about + 15% for Burland’s procedure. The differences 
_ between the measured and predicted values are generally within the engecies 
_ uncertainty for each of the methods. The variability about a mean ‘“‘predicted- -to- 
" measured” ’ ratio is less for Burland’s method and greatest for the API method. 
‘The Burland and \ methods tended to overpredict the capacities. The predicted | 
capacities for the Burland method ranged from about 109% to 123% of the a 
_ measured capacities. The predicted capacities for the \ method ranged from 
“of the in to 115% of the measured capacities; this may be due to an overestimation 
of the in situ overburden stress. The results from th the API method 


| 
J Effect of Loading Rate.—To protect the test piles and loading equipment during | 
the fast tests, a relief valve was installed in the hydraulic system to prevent = 
P a _ overloads. Prior to testing, the relief valve was set to a pressure considerably 


above the maximum pressure expected, during testing. The made 
on a 2-in. (50-mm) diam pressure gage, corresponded to | 408 kips (1,815 kN) 
for Pile 1, 526 kips (2,340 kN) for Pile 2, and 666 kips (2,965 KN) for Pile | 
4. The maximum loads developed during the fast tests were much larger than 
those for the slower pile tests. Maximum loads obtained in the fast tests 
_ corresponded to about 93% of the relief-valve settings. . The difference in load 
between the | > Maximum obtained loads and the loads ‘set on on the | pressure relief on | 
valve may be attributable to a premature opening of the r relief valve, to inaccuracy oe 
in setting the valve, ortoacombinationofthetwo, 
_ Loading rate was defined as the difference between yield load and load at 
_ the beginning of the test divided by the elapsed time between initiation of cin 
mi and reaching the yield load. Fig. 10 shows the ultimate capacity versus s loading — 
_ rate for Pile 2. Although the scatter is large, a definite trend exists. A comparison 
of the fast test results with the first compression test of the second series — 


demonstrates a significant increase in capacity with loading rate. The capacities 


as 
AIG. 10.—Influence of Loading Rate | on Capacity (1 kip = - 4.45 kN; 1 kip/min = 
from the fast tests assume that the measured ultimate rit was a measure > 
of the pile capacity and not a measure of the load associated with an opening 
of the relief valve. If, however, the opening of the relief valve affected the — 


maximum load applied to the piles, the intorpectations on and conclusions | derived 


The results showed that ‘the ‘ultimate pile capacity may be increased by 40% 
to 75% if the loading rate is increased by about three orders of magnitude. 
Even if the measured ‘‘slow’’ compressive capacity, however, is increased by 
_ the full computed end bearing due to the uncertainty in the development of 
- the end bearing, the implied increase in ultimate capacity is not affected a 
‘ significantly. This 40%-75% range between results for the different piles may 7 


_ bed due to differences in the soil properties, differences between 1 the Telative | 
percent of the. load carried in end bearing compared to shaft friction, loading — 
history, or some other causes. Pile 4, e.g., showed the greatest increase in | 
rate. In this case the > fast tests we were re performed ‘Prior 
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after the variable static tests. The less severe loading history for Pile 4 may 
have contributed to the greater in increase in capacity with loading rate. fol 
ed Analyses with simple lumped p: parameter models demonstrated that ms increase 


= at different loading rates showed an increase in undrained shear strength | 
of about 50% for a these- order magnitude increase in loading rate. Application - 
& these results to the: API pile capacity prediction method 1) Suggests | a 


_ loading ra rate. Thus, | the pile test results were in n agreement with the laboratory — 
_ measured increase in undrained shear strength for a corresponding increase 


in rete cr leading baa 


The variable static tests did not reveal any significant increase in —— 


— 


‘the slow tests. The variable static tests, however, were only unloaded to 46%-70% ; 
_ of the static capacity prior to each test. _ These residual seating loads, which 
_ were large compared to the small seating loads of 20% or less for the fast 
tests, probably mobilized most of the soil strength. Thus, the higher loading 
ae. for application of a small part of the total capacity was not able to generate 
‘greater | soil resistance ¢ and greater pile capacity. The variable Static tests were 


ro on Pile 4. The. very ‘small setup | time for n most “of ‘he\vacitble Static tests 
= also explain the lack of any correlation between capacity and loading rate 
_ Effect of Thme.- —An i increase in the setup time should result in 
ieee response of the ‘soil near the pile. ‘Setup time is s defined | as s the ‘elapsed 
- time after driving or the elapsed time since the previous load test, wall 
" The measured pile capacity was plotted against setup time to evaluate the ; 
: ‘{jaftnenes of setup time on pile capacity. The rate of loading has a significantly _ 
- greater influence on the capacity than elapsed time for the conditions of these 
: ” tests. Therefore, only results from the slower rate e of loading tests were included wi 
in the evaluation, and the results of the first test of each series were nesiocted - 
_ for this comparison for reasons discussed earlier. The data showed an —_ 
q in capacity with elapsed time of about 5% per log cycle of time. The impact | 
of variables other than setup time for even t the slower tests is ‘greater than 
that of setup time precluding a any quantitative conclusions on the changes in 4 - 
Effect of Cyclic Loading.—For the applied loading history, the cyclic loading = 
appeared to have more effect on the load-displacement response than on capacity. 
a Consequently, the effect of cyclic loading on pile capacity will be introduced — : 
in the next section on load- ad-displacement response of the 
_ The load-displacement response was evaluated in terms of the secant spring 
constant at yield load for all of the tests and comparisons of the load-displacement 
curves for selected tests. The measured values ranged from about 500 kips/in. 
to 3500 6 to 613 kN/mm). For any tion 


he variahie while for the other nilec the tac : vere Oormed 

| 
= 


in the : stiffness fie the different tests was large. The potential ae oe 
in the measured stiffness in the first series was about 16% to 73% due to i 
: potential errors in the displacements. Modifications were made to the vouamentl 
tation for the second t test series, and these potential e errors W were only 2 2%. é 
Spring constants were predicted with a computer program based on laste 
solutions but with slippage between the soil and pile once the limiting shear — 
; : stress on the soil-pile interface was reached (9). The computed values ranged — i 
from 540 ng ip 6 kN/mm) for Pile 1 to 1,120 kips/in. (196.2 kN/mm) 
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: FIG. 12.—Influence of Loading Rate on Pile-Soil Stiffness (1 kip /in. = = 0.175 am pina; 
Tkip/min = 4.45 kN/min) 10: gaivinh 
Effect: of ‘Loading Rate.—The load- -displacement curves for Test Pile 2 of - 
the fast tests are compared in Fig. 11 with the response of the first compression — 


test of | the second series. These comparisons illustrate a a higher pile-soil stiffness 
plotted against loading rate in Fig. 12 to further evaluate the pile-soil 


stiffness. The degree of scatter is large, but some generalized trends can be — al 


| The pile-soil stiffness of the variable Static tests rally exceeded the value = 
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for the - tests even though the rate of loading 1g was =e the highest. As” 


‘mentioned previously, the variable static tests were only unloaded to about 
55% of their static capacity. The values of pile- -soil stiffness measured in this 
ease are for reloading between about 55% and 100% of the capacity. The reloading _ 
Ww response for a partially unloaded pile is expected to be stiffer than that for 
_ A representative increase in pile-soil stiffness per decade of loading rate was © 
10% for Piles 1 and 2, and 25% for Pile 4. The results of Pile” 3 exhibited 
considerable scatter, and the range in loading rate was too small to determine 
_ any trend. The increase in pile-soil stiffness with loading rate is less than the _ 
‘increase in capacity, which is consistent with theory. For example, if the = 
modulus and strength that can be mobilized are doubled due to loading rate, — 
_ the computed increase in pile-soil ‘stiffness due to ) the n modulus ‘increase was 
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FIG. 13.—Load-Displacement Response under Cyclic Loading ai in. = 25.4 mm; 1 


; # Effect of Elapsed Time. —The rate rate of loading and the m magnitude of the seating 
load have a significant influence on the stiffness. Therefore, only the slow | 
tests were used to evaluate the relation between pile-soil stiffness and elapsed 
time after driving or the previous test. The data showed that the effect of — 
setup time on the pile-soil stiffness was negligible. 
_ Effect of Cyclic Loading.—The load-displacement response of the four piles 
- subjected to one-way cyclic loading are compared in Fig. 13 with the static 
load-displacement curves of the last compression test prior to cyclic loading. — 
The cumulative permanent displacement is shows in Fig. 14 against the ratio — 
of maximum cyclic load to static capacity. Results from the Empire cyclic | 
tests are in general agreement with data reported by Broms (3). The results — 
_ demonstrated for the given load history that piles can support coyorl loads —_ . 


0 00015 


J ‘The: were tend” controlled. ‘When the load reached the p predetermined 

j level, the direction of load was reversed. At higher loads, the rate of deformation _ 
had to be increased to maintain the period. Displacements began to accumulate 

= the maximum cyclic load reached about 80% of the static capacity, and 
these displacements increased appreciably when the maximum cyclic load 

| pace 110% of the static capacity. At very high loads Gis was very 


atin 


st @ LOAD INCREMENT EXCEPT AS NOTED. 
(2) EACH CYCLIC LOAD INCREMENT IS ABOUT 40% OF 
STATIC CAPACITY 


(3) STATIC_ULTIMATE CAPACITIES ARE BASED ON THE 


AIG. 14- —Influence of Leading on Permanent Displacement (1 in. = = 25. 4 mm) 


PEAK CYCLIC LOAD 
STATIC CAPACITY (TEST NO. 2)- 


pene FIG. 15.—Pile-Soil Stiffness under Cyclic Loading (1 kip/in. = 0.175 kN/mm) _ 7 
‘The: pile- -soil stiffness during cyclic loading» was determined for several 
= -displacement hysteretic loops. The measured stiffness, which was defined 
by an average slope of the loops, is plotted in Fig. 15 against the ware 
of load cycles and in terms of the ratio of the peak cyclic load to static capacity. 
The stiffness decreases with increases in the number of cycles, the cyclic load, 
_ Data are available at Bridges 2 and + for the ‘cyclic tests on Piles 3 and 
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4. For Pile "data at Bridge 4. ‘These data were used to evaluate > 

the redistribution of load in the pile during cyclic loading. The load distribution ng 
was initially transferred to deeper pile sections with each successive load cycle 

of a given amplitude, but after several cycles this trend was reversed and the — 

upper portion of the pile began to carry more of the load again. CA a 
_ _The fast tests performed after the cyclic tests gave capacities that exceeded | 

the static capacity. | Thus, the capacity after cyclic loading did not appear to 
be affected significantly by the type of cyclic loading used in this study. —— | 
P-- pile capacity increased with an increase in loading rate, but the effect _ 
appeared to decrease as the of the residual seating load 


"capacity, the increase in n capacity, with an increase in loading rate was 
The cyclic tests, however, had seating loads of up to 90% but were able to | 
_ develop loads as high as 130% of the static capacity and would probably have 
gone higher. The loading rates pila cyclic tests, however, exceeded those 
of most of the variable statictests.§ | 

The loading rate for the maximum loading conditions of most offshore =. 

4 is controlled by ocean wave or earthquake loadings. These loading rates are 

- two or four orders of magnitude greater than the rates for which static capacities 
are based. Consideration of loading rate for offshore piles may have important 
economic implications. As other studies have shown that cyclic loading — 
degrade capacity, the total effects of the loading history must be evaluated 
together with loading rate. _ For the one-way cyclic: loading history used in this 
study, the combined effect of cyclic loading « and rate e of loading ‘resulted in 


Comments 


tests demonstrate the influences that the combined but complex 
"interaction of loading rate, number of load cycles, magnitude of cyclic load, 
loading history, and magnitude of sustained load have on pile aaiiinaion. 
_ Each of these factors influences pile performance in a complex manner, and — 
_ their individual and combined effects can be significant. In view of the complex 
- loading history on these piles, conclusions are more qualitative than quantitative. * 
_ Nevertheless, the results have added to our understanding of pile e performance, | 
especially under cyclic and fast rates of loading. 
_ This research has provided additional insight into performance of short piles | 
' ‘embedded in strong but slightly underconsolidated soils at depth. Results from — 
_ these tests on piles only 40 ft and 50 ft long cannot be used to test conclusively — é 
3 the applicability | of the predictive methods to piles of 200 ft to 300 ft or more 
a in length that are used to support offshore structures. Pile test data on on long 
piles are still needed, as the information generated by this study does not provide 
a complete answer to the performance of long, large-diameter piles. ae we 
tests were funded by 14 oil companies: _Aminoil, British Petroleum, 
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Connell, and Gordon E. Strickland. Soil borings and geotechnical engineering © 
_ were performed by McClelland Engineers, Inc. Field construction was under 
_ the joint control of Chevron and TERA, Inc. Pile instrumentation, field tests, # 
and data were the responsibility of TERA, Inc. 
for Planning, Designing, and ‘Constructing Fixed “Offshore 
Platforms, American Petroleum Institute, API RP 2A, 10thed., 1979. 


_ 2. Bergemann, H., ‘‘Alternating Loading and Pulling Tests on Steel I-beam Piles,”’ 
Proceedings, 8th International Conference on Soil Mechanics and Foundation Engineer- 
ing, Vol. 2.1, 1973, pp. 13-17. 
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_ Pipe Piles in Clay,”’ Proceedings, 11th Offshore Technology Conference, Vol. 2, 1979, q 
6. Doyle, E. H., McClelland, B., and Ferguson, G. H. , “Wire- Line Vane Probe for 
Penetration Measurements of Ocean Sediment ” Proceedings, 3rd 
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_ VIBRATIONS OF EMBANKMENT DAMS 


a _ The types of serious damage that are most frequently o observedinembankment  _f 
_ dams after strong earthquakes can be grouped into the following four categories a 
4 in order of increasing frequency of occurrence (9); 
ae 7 Longitudinal cracks that develop especially near the crest due either to 
tensile stresses caused by the lateral vibrations of the dam or to ‘differential — 
settlements caused by different degrees of dynamic compaction of various dam 
zones (core, shell) or foundation soils. 


* _ 2. Slides of soil masses from the upstream or downstream s slope of the dam 
tr riggered by lateral vibrations that induce along sliding surfaces shear stresses 
exceed, on the average, the shear strength of the soil. 


either due to tensile ‘stresses created by longitudinal vibeatiitis | of the dam or 
due to differential settlements of various sections of the dam in the longitudinal - 


_ 4, Large displacements or - bearing ciplicny failure due to liquefaction of © 
cohesionless soils in the dam or the foundation. han 


Transverse cracking due to longitudinal vibrations is neither as serious ; damage : a 
as are the large permanent deformations caused by liquefaction or sliding, nor x 


_ as frequent a phenomenon as is the development of longitudinal cracks. ‘This 


_ perhaps explains why so little attention has been directed towards understanding — 


the behavior of earth dams when excited in a direction parallel to their axis. — 7 
In contrast, _ extensive experimental and tl theoretical research has been related 
to lateral vibrations of earth dams or embankments (1,2,3,4,5,7,9). 


4 ‘Asst. Prof. of Civ. _Engrg., Case Western Reserve Univ., Cleveland, Ohio 44106. . 
_ Note.—Discussion open until June 1, 1981. To extend the closing date one month, 
4 written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on December 18 a P 
; 1979. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 


of the American Society of Civil Engineers, © ASCE, Vol. Wi, GTI, 
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However, ‘transverse ‘cracks can be p potentially very. since they 
iat can induce uncontrolled seepage through the dam and thereby create piping J 
failure(s). Such cracks are highly probable to occur especially at the contact 
surfaces of an embankment dam with its steep-sloping rock abutments, because _ 
‘4 of the different dynamic response of the two media, and because no tensile _ 
‘om _ forces can develop between them. The 1/16-in. wide (1.6 mm) lateral crack 
x, which occurred at the east abutment of the Santa Felicia Dam during the 
Fernando Earthquake of February 9, 1971, as reported ples a oe and 
Scott (1) seems to have been of thismature, 
The contact areas between soil and outlet works such side Spillways, 
2a sluiceways, and buried pipes are also vulnerable to tensile longitudinal stresses. 
eZ Furthermore, cracking and ruptures of pipes observed after earthquakes have 


been attributed to the incompatibility of their displacements with the larger a 


- occur within the body of the dam and can easily Temain undetected they may © 
cause additional problems in the post-earthquake life of the structure. Conse- _ 
quently, experimental and theoretical research on the behavior of embankment ry: 
dams during earthquake- -induced longitudinal vibrations is needed to to 


homogeneous ‘triangular’ prism, bounded in n the longitudinal direction by t two 
: SS planes (rectangular canyon). Numerical results are presented for the 
= frequencies, modal displacement, modal strain shapes, and modal parti- 
__-Gipation factors. It is shown that in ne ge long dams shear deformations 


1% in 1 order to evaluate the ability of the presented theory to explale with reasonable 
observed field behavior. Natural frequencies, participation 
a and mode shapes computed with the theory for the Santa Felicia DaminCalifornia | 
and the Kisenyama Dam in Japan are examined in the light of observed 
predominant frequencies and peak Tesponse amplitudes of motions recorded * 
_ during actual earthquakes (1,9). Finally the method is extended to approximately 
account for a realistic of shear modulus at various below 
8 crest of the dam. dig 


hag 


ANALYSIS FOR For Homocencous | Canron 
wg Free Vibrations.—Let x, y, and z be the ecthueadeials evestianbes of a any va 
” in a symmetrical earth dam built in a rectangular canyon, as shown in Fig. | 
- I(a). To derive the governing equation of meme | the simplifying assumption 
is made that normal and shear stresses, o, and t,,, are independent of y and, 
bs ani therefore, that o, and t,, are uniformly distributed over infinitesimal areas b,-dz 
and b,: -dx respectively, as shown in Fig.  1(6). assumption is reminiscent 
cd of the hypothesis of uniform shear stresses T,, over areas b,-dx ‘that is the 
basis of the for (3,4,7,9). It seems, 


at 
| 
| 
wa design procedures, and invent rehabilitation techniques to ensure the safety or 
Of new or existing earth and rockfill dams during and after earthquakes. 
1 a _ Asa first step in the aforementioned direction, this paper presents a theoretical sg 
tag » 
analysis of free and forced vibrations of embankment dams. The method accounts ff 


faces of the dam, in contrast with the shear- beam assumption. ag 


‘As a direct consequence of the foregoing assumption, the horizontal 


ments are independent of y, i.e., u z;t). The differential 


x dx 


> 
dx “3 


derived the dynamic of an body of volume b, ~ 
in which G and E = and Young’ s moduli Py the soil, herein taken 
as constants the dam. The three terms in Eq. | express 


a 
TTT] 
| and (8) Longitudinal and Cros 
1.—(a) Dam Geometry 


Solutions: of ‘Eq. physically only if they satisfy the 
ern 


La 
that e express the requirements of zero relative displacement between dam = : 


dam and 
u= = ¥(x) exp 


and call k = w/C, in which the oars v 
‘operations Eq. | becomes 


Each of the terms in Eq. 


l is a constant, a’ 


of z z only and the second term a a function of x only. 
+ +——V = 


general of Eqs. 8 and 9 that ‘satisfy the boundary Eqs. and 


= Oey = B” 


13) 
in which = integration constants; Function of the 


elocity is C = [G/p]'’*. After some § 
4 
@ 
| 


= 8.6530, etc. 
The natural tre requencies for longitudinal vibrations a are ites by > 
while ‘displacement table. normalized to a u are 
Bq. 240 states that, for a 


(+ <) ‘sin 


’ 


the normal and shear strains 


bd 


dane of vertically propagating S-waves polarized in the x direction, the _ 


- longitudinal motion of the dam relative to the surrounding canyon is described = 


gig 
in in which i,.(t) = = the ground acceleration (in ‘the x direction), and by Egs. 
2, 3, 4, and 5 (boundary — The — obtained by modal 
superposition (8): 
ton Bis. | 


to, the modal participation factor e tual the n, r mode is 


Ales 
ish 
“while ‘modal response D,,,(t) is obtained from the Duhamel ‘Gon 


wet 

~ ax Ox 


in which ne = damping ratio ratio in the n, r mode. 
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“Fig. demonstrates the relationship between natural frequencies, of a 


dam in a canyon and the re geometric parameter, 


© ab adi to. lanibooly 
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«FIG. 2 —Dependence of Natural Frequencies on Length-to-Height Ratio 
“y of Poisson’ s ratio, v = 0.30 and v = 0.5. One can clearly distinguish two regions | 
in the frequency-ratio spectrum, corresponding to ** narrow’”’ or ‘‘wide’” 
_ Tespectively. When L/H is less than about 1. 5 (‘‘narrow”’ canyon) w,,H/C 
AS is inversely proportional to L/H, whereas when L/H is greater than about 
2 (“wide” canyon) w,, H/Ci is of L/H. The following 


: 4 of the fundamental frequency-ratio, w,,H/C, as L/H increases, fortwo values [im 
ie 
fa 


ongitudinal of “narrow or “long” dams: od? 
vlog 
sh om ai swords 


Eq. states for a dam, T,, is independent the eight, 

and is only a function of the length, ond the S-wave velocity, 
esult can be compared with ‘the fundamental period in axial vibrations ofa a 
uniform rod fixed at both ends and having a Poisson’s ratio v= = 0. val ee ed 


: The similarity of Eq. 24a and Eq. 25a suggests that axial deformations are _ 
more significant than shear deformations in narrow dams. a 


a ‘The reverse appears to be true in long dams; T,,, depends only on the height, 
a of the dam and Eq. 2 is reminiscent of the fundamental period of a two 


wedge-shaped beam due to lateral shear vibrations (3,4,7): 


shear deformatio ore than axial deformati 


‘which also 4 also demonstrates that as the canyon becomes wider (H/L decreases) 


_ the importance of shear deformations increases. 7 

= Notice also in Fig. 2(a) and in Eqs. 24 that Poisson’s ratio, v, , has 2 a a practically =. 

— negligible effect on the fundamental frequency. Also negligible is the influence — 

_ of v on the higher frequencies. Only results for a typical value of v = 0. 40 

are thus shown hereafter, except when otherwise | 

-_, Fig. 2(6) portrays the variation with L/H of the ratio on te of six higher 3 
frequencies (n, r = 1, 2, 3, 4) over the fundamental frequency. Two :comienons = : 
can be drawn: (1) Successive natural frequencies are very close to each other, 
e.g., a dam with L/H = 2 has w,, = w,, = 1.60 = 1.40 

- = 1.06 w,,, etc., therefore, in general, no reliable estimate of response parameters Lal 
can be obtained through empirical re of the maximum modal 


n spectrum (8). Instead, 4 


i! 
> 
(255) 
ons in j 
_ deformations is by considering the peak-strain ratio == 
(SR), = 3.10 A | 
_ | of time histories of modal responses (as 


0) 


the ratio w,,/ is an or sanction of L/H depending on 
_ whether n is larger or smaller than r. When n = r, the ratio [not shown in @ 
“Fig. 2(6)] has an approximately constant value, i.e., independent of L/H (e.g., 
b Finally, F Fig. 3 shows > re the eee, 20 of an earth dam in the 


3.—Distortion of Dam Vibrating in Two Natural Modes (perspective) 
MODAL DISPLACEMENTS 


> 


| 


i 
MODAL "STRAINS 


“FIG. —Modal and Modal Strain of First Two 


and First Two Antisymmetric Modes 


first two symmetrical modes of longitudinal vibrations [(1, 1) and (2,1)] and 
_ Fig. 4 displays the x and z components of modal displacement shapes with 
_ the corresponding axial and shear components of modal strain shapes and the 
mode participation factors. It. _can be seen that maximum axial strains (which 


are of major concern with longitudinal vibrations) develop at the crest level — 


of the abutment dam interface. etiienmes contributions to these strains come 


a 
| 
al) 
| 


7 aa = 0.75 H or 0.34 H etc. depending on the mode of vibration (n ‘= 


sia 
Abdel- -Ghaffar and Scott have reported two “complete” 


tke 
god 
| between pre pene 


UPSTREAM 


¥ the San Fernando Restheuabe: of February 9, 1971 (M = 6.5, R = 
the ‘Southern of (M = 5, R= 


EMBANKMENTDAMS 

— from both symmetric (r = 1) and antisymmetric (r = 2) modes of | 
7 
| 

km) 
4 


_characteristics are shown in Fi ig. 5 Secccad from Ref. 1). The dam, on ll 
; 65 km northwest of Los Angeles, is approx 83 m high above its lowest foundation, — 

61 m above the original stream bed, and 390 m long at the crest. It has a _ 

central impervious core and shells Upstream and downstream. All 


TABLE 1 between Computed and Observed Natural Frequencies and 


Participations: Santa Felicia Dam, California 


San Fernando” . California Co 

Earthquake: rin | Earthquake: -4 = 278.3 m, H = m, 
_February 9,1971 |  April4-8,1976 | 230,» = 0.45) q 


Mode 
in hertz 


| 
340 


Oe YUWN We 


’ Following the two earthquakes, comprehensive full-scale dynamic tests were 

out, inc! tuding m mechanical, ambient, and 1 hydrodynamically-it -induced vibra 

4 
 low-strain ‘values of key soil parameters such as shear-wave velocity (or shear 

pene and Poisson’s ratio (1). Wave velocity estimates were also obtained q 
from observed resonant frequencies of lateral vibrations, utilizing existing 

shear-beam theories (3,7). It was suggested that the average velocity of the | 


in Mertz in Nertz | A/A 
1350 | 1.00 | 1.00 1411 1.00) 
| 038 2893 | 
| 
322 | | 3.140 0332 | 2, 
= | — } 
4.2000 (0.79 | 3.975 | 0.166 | 2, 
442 | 036 | 459 | 0.44 4,456 0.531 | 3, be 
=: 4501 | 0.133 | 28 
ji 0.83 9.838 | 0.153 | 68 
= peak amplitude of the crest /abutment amplification 


Santa Felicia dam lies between approx 220 m/ sec and 270 m/ sec, while 
Poisson’s ratio equal to 0.45 seems reasonable. % 

_ An “‘equivalent’’ dam in a rectangular canyon shown in Fig. 1(a) is chosen — 

to approximate the actual dam geometry. Its dimensions are: L = 278.3 m 

oF and H = 72 m. Shear-wave velocity is taken equal to 230 m/s, i.e., near the _ 

x lower bound of the suggested range, as is ; appropriate for “relatively strong — ; 


earthquake motions that are analyzed here. movi diqsb dtiv 
The ratios of Fourier i.c., computed from 


“the dam i in longitudinal vibrations. They also offer a a qualitative picture of pq 
degree to which various natural modes participate in the shaking. Table _ 
summarizes this information (obtained from Ref. 1) in the form of natural 
Ps frequencies, f, and ratios of peak amplitudes, A/A,,. Also shown in this table 
B are the natural frequencies and mode participation factors computed wi with the 
theory presented herein (Eq. 16 « and Eq. 22) using - the aforementioned 
for the material properties and geometry of the dam. Se ee 
The agreement between predicted and observed natural frequencies is a 
all purposes satisfactory although some of the predicted higher frequencies are _ 
not observed in the recorded motions. Notice, however, the relatively small © 
participation factors of these 
_ Direct comparison between ratios of peak amplitudes (A,,/A,,) of an o engl 
fication spectrum and ratios of participation factors ([,,/T,,) cannot be made 
_ Because as it can be seen from Eq. MRA dai “perl 
and, unless i(t) |max» the aforementioned two ratios will ibe 


different. Nevertheless, one can qualitatively state that the relatively large 4 
| Os participation factors predicted by this theory for some of the higher = ' 
_ of vibration are in accord with the large peak amplitudes « of the amplification — > 
4 In conclusion, the presented simple ‘theory of longitudinal vibrations © 2 
_ ‘homogencous earth dams in rectangular canyons explains with reasonable 7 
7 accuracy the observed behavior (natural frequencies, participation factors) of 


(9) reports the motions recorded during a 1969 Earthquake by a a 
BP ero of seismometers that had been installed at the crest, slope, interior, — 
and abutments of the Kisenyama rockfill dam, in Japan. The dam is 95 m 
high, 255 m long at the crest, and is founded on rock. Fig. 6(a) shows a longitudinal = 
and vertical cross section o of the _ dam, while Fig. 6(b) shows the parallel-to-the- 


dam-axis motions recorded at the ground (point 1), at t the center of the crest , 

(point 2), and at 25 m below the crest in the clay-core (point 3). The earthquakes © 

occured at the northwestern part of the ioum Prefecture on ry 9, 1969, 


| 
recorded at the middle of the crest and the right abutment during the <a 
| 


— 


modes of the dam are excited; the motions (at point 2 and point 3) are very. 
nearly sinusoidal with an average period T = 0.42 sec. Also apparent is the | 
huge amplification of the motion at the crest (point 2), whereas at a relatively _ 
shallow depth of approximately (1/3) H from the crest (point 3), the motion 
is only moderately higher than the ground motion (point 1). The distribution 
with depth from the crest of the peak recorded accelerations is depicted in 


> 


rectangular canyon 


0.08 


location 


a FIG. 6. aiteabitte Rockfill Ted (a) = (b) Accelerograms i in Dam and on 


2 and (c) Variation of from youre 


sec; and 2) large amplification of the motion near the crest. 
Estimation of S-Wave Velocity.— —On the basis of the observed fundamental 
period T = 0.42 sec, the | theory of longitudinal vibrations of embankment dams 
"presented i in this paper can establish an ‘‘effective’’ value of the average S- wave 
velocity, consistent with the experienced level of shear deformations during 
_ this particular earthquake. This velocity can then be utilized to explain the 


the records reveals that essentaliy only ine tundamental and a few hisher 
| 
location of seismometers P 


EMBANKMENT DAMS 


“actual geometry, as sketched in Fig. 6(a). Its dimensions are L = 168 m ae 
as H = 82 m. Poisson’s ratio is taken as v = 0.25. This is believed to be a _ 
representative value for rockfill dams immediately after ccastruction, i.e., before 

the water of the lake has saturated the clay core. Sar 

s Now, substituting T,, = 0.42 sec and the aforementioned geometric cand material 
"properties in the frequency relation (Eq. 16) gives; 


( 


©. 42) 168 / 


or C = 359 m/sec. The theory of lateral shear vibrations of embankment dams 


ina rectangular canyon (3,7) is subsequently used in connection with 1 the observed — 
fundamental period T= a 0.48 sec to check the accuracy of the preceding value ae 
"velocity. Solving for CG, the frequency e.g., = 


or C = 375 m/sec. The | backfigured values of Cc are in fairly good agreement 
with each other (difference ~4%). Moreover, 359 m/sec-375 m/sec is a realistic 7 
_ Tange for the S-wave velocity of a rockfill dam experiencing moderate levels = _ 
5 4 of excitation (~0.10 g peak-crest acceleration i ‘in the lateral direction; and ~0.07 = 
‘The natural periods of the two higher symmetric modes of vibrations are are 
obtained means oe (6) (or through Eq. 16): intern! vit whats 


0b) 


- One can easily identify these two modes with some of the peaks of the recorded 
motions, e.g., the crest accelerogram (point 2) at time ¢ mali sec and in 7 


_ have a negligible contribution to the motion of points 2 and 3 _ because — 
are located close to the central cross section of the dam, i.e., at x= = L/2. 

_ Distribution of Peak Responses along Height of Dam. —To explain the sharp 
increase of the motion near the crest of the dam in comparison with the motion 
at a depth of only 25 m (=0.305 H) from | the the crest, one can use the presented 


0.20 sec. Note, also, that the higher antisymmetric modes [(1,2), (1,3), etc.] : 


| 
| 
| 
| E 
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‘Thea tiene histories at points 2 and 3, ii, (t) and ii, (t), canbe 
_ obtained by means of Eq. 21 after evaluation of the Duhamel integral (Eq. 
23). Unfortunately, however, the ground record u,,(t) is not available in digitized — * 
form and in order to roughly estimate AR the response spectrum approach ~ 
- (8) is used with the following two > simplifying assumptions: (1) Only the first 
z three symmetric ‘natural modes” (n = 1, 2, 3; r= 1) participate in the m ‘motion; 
and (2) the spectral accelerations, S,, that correspond to each of the three 
Z first natural periods are all equal. The first assumption is by no means vs 


as it is evident from the appearance of the two accelerograms and the previous 


Modal Maxima HOMOGENEOUS" THEORY "INHOMOGENEOUS THEORY 


4 i 


rr: s s 4 25 


— 


Rove 
ASSUMPTION: )=s a(T3)) = * = 2° 


» 


w ot 


i 


_ analysis regarding the natural periods of the dam. However, since the oo 


of the response spectrum are unknown, the second assumption can only be 
_ considered as a crude, although reasonable, simplification. Its justification stems 


from the observation (8,9) that usually damped acceleration spectra exhibit a 
more or less constant value in the period range between 0.16 sec and 0.42 
sec, i.e., in the range of the three natural periods of interest. This would especially 


be so for the large value of critical damping ratio (of the order of 5%- —_ 


q 
| “TE | 
q 
For each of the two points (2 and 3) the maximum acceleration due —. 


EMBANKMENT 


0,= max ii, r BOs. L 2, ¥ to ano) 
Ba he Soi 1, est (2,0) od? (3,2) 
4 in which the participation factors, T,,,, and modal shapes, U,,, of the (m,1) 
mode are . obtained from Eqs. 22 and 17, “respectively. ’ The individual modal 
maxima are then combined by three empirical rules to provide an estimate 
a of the overall expected peak accelerations. Table 2 displays the results and 
. explains the expressions of the three rules. The “‘square root of the sum of _ 
am the squares’ (SRSS) is the most popular of these rules (8,9) but it eed 
~s overestimates the total-peak response in this case because of the proximity 7 
2 7 of the natural frequencies of the dam (a phenomenon a already addressed in 
E a previous section). The second rule, suggested by Rosenblueth (8), is perhaps _ 
more appropriate since it accounts for the interaction between modal — 
tions. Finally, the third rule attempts to reflect the ‘greater importance of the i 
fendamental mode, as seen in the two records. 4 
The results of Table 2 demonstrate that the simple theory of longitudinal 
é vibrations, as s presented so far, cannot adequately explain the observed high | 
4 amplification of the seismic motion near the crest. All three rules yield quite 
P similar acceleration ratios not exceeding 1.44, which is only 60% of the ratio 
of recorded — accelerations: 0.68 /0. 29 = 2. 34. Of course, one might argue 
0.16 sec) would have favored the “corresponding higher modes ‘and, ‘thus, led 
to prediction of sharper near-crest amplification. Nevertheless, even if “one 
assumes S,(T,,) = S A(T) = 2:S AT )s the resulting acceleration ratio hae 
not exceed the value of 1.88 (shown also in Table 2). Larger differences in 
the spectral accelerations (than the factor 2 implies) are highly improbable given 4 
_ the relatively large amount of damping the dam must have experienced during 
7 the earthquake. After all, the two records bear no evidence whatsoever of 
a any predominance of the higher modes. Note, furthermore, that lateral vibrations 


exhibited a similarly sharp amplification near the crest (8), as pga 
oe although the corresponding natural periods (0.48 sec, 0.235 sec, 
and 0.158 sec) differ somewhat from the periods of longitudinal vibrations. hee 
i, In conclusion, the simple theory that was previously presented cannot explain 
gs distribution of peak responses along the height of the dam. It is believed 
the soil inhomogeneity, which the theory does not account for, is the — 


Mooet tHat Accounts For Son INHOMOGENEITY 

- Shear modulus even within a uniform mass of soil is not constant but increases _ 
approximately as the square-root of the effective normal octahedral stress: G_ 
Plane-strain finite-element analyses of the static stress distribution 
in typical dam cross sections (Ref. 5) suggest in a reasonably consistent way 


that, with good accuracy, the average a4 across the dam width is an | increasing + 


‘Poth’ 


| 
4 
A 
increases as the two-third power of the 


oenesinteniaiae and rockfill dams in the United States and 
Japan have confirmed this form of variation of soil stiffness in embankment-type _ = 
dams (5,6). The writer (4,5) has recently presented a theory * sateral shear 3 
_ vibrations of earth dams whose modulus varies i in Proportion to z’’*. Successful ce 
- evaluation of this theory in the light of recorded responses of fhe dams 
1] from United States, Japan, and Yugoslavia during earthquake, forced, and ambi- — 
_ ent vibrations (5,6), offers additional indirect evidence of the mp er el 
st is concluded that embankment dams are inhomogeneous. in ‘the nea 
een, It would, therefore, be necessary to modify the presented — 
of longitudinal vibrations in order to account for an S-wave velocity that increases 


= the velocity at the he base of of the dam, i. €., at at A. 
“S Exact analytical solution of the g governing ng differential equation of free vibrations — a 


a constant Poisson’ s ratio. A simple approximate solution can, however , be 
ee irne: by heuristically combining the solution of two extreme cases: (1) The © 


case of a ‘“‘very’’ long dam; and (2) the case of a “‘very’’ narrow dam. As 


explained previously (Eqs. 24 and 26), shear deformations are predominant in 


case whereas axial deformations are re predominant in case de Accordingly, 


‘with G= = (6/7)? > being the average shear m 
exact of Eq. 35 that satisfies the conditions (Eq. 4 


cm 


w 

in which C = (6/7) C,, is the > average S-wave velocity. The: reader is referred 
to Ref. 4 for a detailed derivation of this ‘solution. Direct ‘substitution of Eq. 
37 and Eq. 38 in Eq. 35 is, however, sufficient to convince of its correctness. _ 


a 
q 
affubute 


— heuristic fashion, guided by the form of Eqs. 16 and 17, ‘one can combin : 
Eqs. 37- 39 and Eqs. 38-40 to obtain approximate general formulas for the 


- Themodalpa modal participation factor due to earthquake shaking that consists of vertically : 


pagating waves is also given here: 


the... 


in the two limiting cases of L/H — 0 (narrow dam) | and of L/H — o (long 3 
- dam). It is moreover believed that for intermediate ranges < of the L/H ratio, 
Eqs. 41 and 42 describe the behavior of the dam more accurately than Eqs. = 
16 and 17, since the former capture both the two dimensional character of r) 
_ the response and the inhomogeneous nature of the material properties. » avorest- 

Comparison of Eq. 42 and Eq. 16, portrayed in Fig. 7(a), reveals that 4 

_two theories (hereafter referred to as “‘homogeneous’’ or “‘inhomogeneous”” 

predict ve very similar natural frequencies. The maximum discrepancy 

“tends toe frequencies is only 1. 6%, in the extreme case of L/H = ~; it 


4 
Notice thet its approximate character, the solution i is exact 


tends to zero as L/H decreases. Larger discrepancies are observed with the _ 
higher modes, as the inhomogeneous theory predicts frequencies that are closer 
to the fundamental frequency, in comparison with the higher frequencies" of 


the homogeneous theory a. 16), e.g., for a very long dam Eq. 42 yields 


a however, L/H decreases, the discrepancies become insignificant “nuit 
(e The important difference between the two theories lies i in their mode shapes. 
AS. cows in Fig. the inhomogeneous theory (Eq. 41) leads" toa ‘sharp 


in 
g 
4 
a 
: 
¥ 


| 


of modal not by the 
_ homogeneous theory (Eq. 16). This discrepancy is not decreasing with L/H, 
although its consequence on the shear- Strain /axial- strain ratio is, indeed, dimin- | 

ishing as the dam becomes narrower. 


ono bow AL apa to aa ott yd . wi 


inhomogeneous model > 


, Of homogeneous mode 


ere, 
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7. of: (a) Three’ Natural Frequencies: and (b) Vertical Components | 
of Two Mode Shapes Obtained with (Eqs. 16, and 


: * First, using Eq. 42 instead of Eq. 16, the effective S- wie velocity of the 


— 
4 
| 
1.04 
vant 
owt 


4 
oo The of: the higher symmetric modes are now expected 
to be 6(a)} b ai | Sab 
a 


0.19 sec (466) 

pared with the values of 0.25 sec and 0.16 sec derived from Fig. oes | 

_ Finally, the peak accelerations at points 2 and 3 are reestimated on the basis 

of the same, aforementioned assumptions, i.e., only three "modes participate 
P in the motion; the corresponding three spectral accelerations are the same. The 4 

. “results are displayed in Table 2 for all three rules of combining modal maxima. 
_ The agreement of the predicted peak acceleration ratios with the observed value _ 

3 2.34 clearly demonstrates that inhomogeneity (of the form described by Eqs. 7 


“the motion near the top of the dam. Also shown in 1 this table is s the prediction 
. ‘the peak values and their ratio under the assumption that the spectral 
acceleration, S,, of the second and third higher modes is two times larger than 
“the spectral acceleration of the fundamental mode. The result (AR = 3.57 > 


\. three modes. Thus, it appears that the homogeneous theory cannot even gold 

explain the observed distribution of accelerations along the height of the dam. 

_ The superiority of the inhomogeneous theory, in this respect, is clear. Whether 
Or not the theory will always yield souuits comparing as favorably with reality 
nd _ as in this Single case, remains to to be s ‘seen . The reader is is } cautioned, nevertheless, 


cases. Their effect should, therefore, at least qualitatively be assessed when 

studying the longitudinal as well as the lateral vibrational behavior of earth 


strong, earthquake-induced vibrations, may also play a predominant role in sas 


engineering since they can cause transverse cracking at the abutment dam or 
the dam-outlet works interfaces. Such cracking is quite dangerous as it may 
allow water to flow through the dam and thus lead to piping failure. 
and forced vibrations of embankment dams i ina rectangular 
a method that accounts { for both shear and 

~ dilatational deformations and | models the dam as a linear homogeneous medium. 

_ Presented numerical results demonstrate the effect of the length and height 
of the dam on its natural frequencies and its modal displacement and strain 
_ shapes. The method has been successfully evaluated by comparing the predicted 


natural frequencies for the Santa F eee Dam with the ow predominant 


& _i.e., smaller by only about 0.84% than the value 359 m/s derived previously J 
4 
Longitudinal vibrations of embankment dams are of concern in geotechnical ; 
7 


frequencies of the amplification spectra ra obtained from recorded motions on 


this Dam during twoearthquakes(1), aon 
-- Failure: of the method to explain the observed sharp increase of the motion __ 


‘near ‘the | crest of the Kiseneyama 1 rockfill dam in Japan during a moderate — ' 
(9) “motivated the study of the importance of soil inhomogeneity. 

_ Anapproximate method has thus been developed that considers the shear modulus 
in the dam as increasing with the 2/3 power of depth, a variation that has 
been confirmed from in situ measurement in several dams (1,2,5,6). This simple 
**inhomogeneous’”’ method is negligibly different from the ‘‘homogeneous’”’ theory 
when natural frequencies are compared. It yields, however, mode shapes oo 
“exhibit a sharp attenuation with depth near the crest, and, consequently, 
explains very well the observed behavior of the Kisenyama Dam. 
In conclusion, it appears that additional theoretical and, even more, experimen- 


tal research may be needed in order to develop reliable methods of design 
embankment dams to resist longitudinal | 
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Cyciic AXIAL RESPONSE OF SINGLE PLE 


» By Harry G. Poulos,’ M. | 


Intropuction int festace of the ts can 


was { from tests on piles subjected to ‘ “one-way” ’ cyclic loading (that is, 
between zero minimum load and a specified maximum load), and these tests — 
indicated that there was no significant effect on pile capacity as a result of 
‘the previous cyclic loading (5). d fiction will frags 
a series of model tests on 
piles in - remolded clay which indicated that considerable loss of skin —— 
_could occur if a pile was subjected to “‘two-way’ cyclic loading (loading alternating 
od between tension and compression with a zero or small mean value), if oe 
load magnitude approached the static pile capacity. Their tests also confirmed — 
previous finding that ‘‘one-way”” cyclic loading little 


indicated a trend toward decreasing load capacity and ‘stiffness with increasing 
number of cycles and cyclic load level. 
Because of the possibility of ‘‘two-way”’ cyclic loading on pile 
a supporting offshore structures, it is essential to develop a means for predicting 

the circumstances under which cyclic loading will affect pile performance. Sangrey 

_ (17) has suggested a method of application of critical state concepts to calculate = 
the ultimate axial capacity of  cyclically-loaded piles, but this approach does 
not consider deformations. Matlock and Foo (12) have described an analysis 
7 4g which a pile in a relatively simple hysteretic and degrading soil model is 
considered; soil degradation is assumed to occur if a full reversal of yielding — 
“4 in both directions occurs. Boulen, et al. (4) have described a ‘finite element — 


analysis in which a hyperbolic relationship between strain and number of cycles 
_ employed to analyze a pile in sand. Poulos (13) has outlined an effective 
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oe develop excess pore pressures in the soil adjacent to the pile with a consequent a, 

reduction in soil modulus and skin resistance. However, some of the : assumptions 7 
-:. which the effective stress analysis was based are of limited generality, _ its 
and the results presented may not be applicable to a wide range of problems. a 
_ Moreover, recent tests by Grosch and Reese (6) suggest that cyclic Copatenen 

_ of skin friction arises primarily from destruction of interparticle bonds, and 

particle realinement, rather than the development of pore pressures. 
A further evolution of the previous analysis (13) is described herein, allowing, 
for the consideration of cyclic degradation of the soil in terms of total oad 


The major objectives are to answertwo questions; 


a 


sof load increases? 


of allowing for cyclic degradation of the soil are described. A series of die > 
: for ¢ a typical offshore pile is then presented and the | effects on n cyclic | axial ~ 


_ _ rps scovens of various pile and soil parameters are examined. Finally, the the 


a The static loading analysis has been described ~ Poulos (3) and “will ‘be 
only briefly reviewed here. It is a simplified form of boundary element analysis 

in which the pile is discretized into a series of elements and the soil is considered 


to be a linearly elastic continuum. To allow for piles of nonuniform section, 


_ three types of elements can be used: cylindrical shaft elements, annular base 


| 
| | 
i 
capacity? 
ge as the number of cycles 
| a 


and elements at diecontiqnities in shaft diameter (see Fig. 
1). By considerations of equilibrium and compatibility between the vertical 
movements of pile and soil at each element, relationships may be derived between | 
gy pile-soil interaction stresses, the pile v vertical displacement, and the increment — 


of axial load. These relationships may be solved for the interaction stress 
canbe ob and the incremental pile deflections, from which the over-all values — 


can be obtained by addition to the existing values. ts” 

_ Allowance for pile-soil slip or yield of an element may be made by specifying Ms 
an upper limit to the value of the interaction stress at each element, and carrying — 
- an iterative analysis. ‘This allows a | Static load-settlement curve t to failure — 
to be computed, and also allows for unloading of a pile to be analyzed. if xt 
-pile- -soil slip has occurred, unloading to zero load will result in residual stresses - 
being developed in the pile (in addition to those which may be present due sae 

_ An important feature of the foregoing analysis is that it can take into account, 

approximately, nonhomogeneity of the soil along and beneath the pile ‘shaft, 

as described in Ref. 14. The ability to consider soil nonhomogeneity i in a convenient 
and economical manner is essential when considering a cyclically loaded pile, 4 
.- “‘degradation’’ of soil properties will occur at different rates along the pile, ae 


Ve 


= on the local stress level. Even in an initially homogeneous soil, a 
/— nonuniform distribution of soil modulus and skin friction will result from cyclic 


_ loading b becat of the nonuniformity of stress distribution along the ‘pile. — fre 
EXTENSION oF Anatysis To Cycuc Loapine 
The Static analysis can be to in ways: 


i? By carrying out an incremental analysis: in which the entire cyclic load 
_ history of the pile is applied, with the soil parameters being progressively adjusted 
at each increment to reflect the effects of the cyclic loading. , OF 
& By carrying out an analysis to ascertain the response of a pile after a 
number of cycles N of uniform magnitude (maximum P,,,., minimum P,,,,) have 


been applied to the pile, with the soil parameters being adjusted - to reflect — 
the effects of cyclic loading at the end of the load sequence. 


The first approach, | although ‘more satisfactory, would involve ‘enormous 
i. computational effort if the load sequence involved a large number of cycles; 

consequently, the second approach has been adopted. pa 
The detailed analysis procedure is as follows: to 
4 1. First estimates of soil modulus . E,, ultimate skin friction t,, and ultimate 

_ base resistance p,, are chosen for each element (e.g., the values for static 

43 2. The Pile is analyzed for the 1 maximum load and and the dist distribution of 

| 

a 3. The pile i is s similarly analyzed for the minimum im load ee eee - 
The cyclic shear stress and cyclic displacement p. are ‘determined 
the two fi oregoing analyses by subtracting the minimum vetoes from the maximum Ps ; 


shear 1 stress, the following. determinations are made: Cyclic displacement < ‘or ¥ : 
_ cyclic shear strain at each element (depending on the type of degradation cad 
_ being performed); the change in soil and ultimate skin friction 


the estimated values, and if the difference i is a specified tolerance, 
new values are selected and steps 2 to 6 repeated until the desired dégree a 
7. The cyclic deflection, the mean aateeniie and the values of ultimate skin — 
friction and base resistance are thus obtained, from which the available ultimate — 
capacity after cycling may be calculated. 


to the value for static loading. Three of degradation 
factor are relevant here: 
‘Ultimate skin friction degradation factor D,,thatis © galbeo! 


‘in which 


b 
_ in which p,,.. 


“in which | E, = = modulus at after cyclic loading; and E, = = soil 
_ modulus at element, for static loading. Normally, the degradation factors will 
“be less than or equal to 1.0, although values larger than unity may apply in 
¢ some cohesionless soils, or if the effects of rapid rates of loading are taken 
The major problem in analyzing the cyc cyclic re: response of a a pile is to determine 
‘how these degradation factors vary with increasing stress or strain level and - 
amber of Some possible approaches are described subsequently, 


= 
4 
; _ To simplify the computations and give a better insight into the progress Of 
_ cyclic degradation, the concept of degradation factors has been introduced (13), 
| | 
to cyche loading, and t,, = ultimate skin Iriction for stauc loading. 
2. Ultimate base resistance degradation factor D,, that is 
q 
ot 
timate 
q 


Mer ETHODS OF STIMATING DEGRADATION now FACTORS 

a The ideal data to have available would be the variation of the three degradation — 

_ factors D,, D,, and D, with cyclic interface shear stress level, number of | 
cycles, poy cycle frequency. Unfortunately, no such data are currently available — 

2 it is therefore necessary to resort to approximate means of of estimating the 


factors. Two bre broad ad approaches could be adopted: pevemeter 


total stress aan in which cyclic loading is assumed to occur under | 


elas conditions, and no attempt is made to distinguish between the various ; ; 
factors contributing to cyclic degradation. This approach is straightforward, 
but has the limitation that no account can be taken of time effects, as conditions 


assumed to remain undrained. flog the effec of 


a (3). By determining the effect of a change in excess pore pressure (and hence 
‘effective stress) on 7, > Po,» and E,, the corresponding degradation factors may 
fs obtained. Such an approach has the advantage that time effects resulting | 


from pore pressure dissipation during cycling may be incorporated. However, — 
= ly 
no allowance is made for other factors contributing to cyclic degradation, €.g., 
; structural breakdown of the soil particles, creep of the soil skeleton, and as 2 
previously mentioned, such factors: may be the cause of the 
4 i total stress approach will be used herein, and it will be assumed that at cyclic 


_ observations of Thiers and Seed (18) and Lee ‘and Focht ( l 1), that the eo 
in undrained shear strength s, is a function of the cyclic shear strain Ye prior 

x _ to failure, normalized with seapect to the static shear strain to failure, y,. Since - 

4 it is reasonable to assume that +, and Pou both depend on s,, it should also 
be possible to Telate the degradation of t, and p,, to y./¥, . However, itis 
well known that full mobilization of skin friction 0 occurs at on wl or or displaceme: Ps 
which are only about 10% to 25% of those associated with either the Ewa 


| 
Gegradauon is dependent On tie cycue Stain level, and Not On permanent 
level. Data to support the validity of this assumption have been presented 

Degradation of Ultimate Skin Friction and Base Resistance.—Because of the — a 
— lack of direct data on the degradation of ultimate skin frictiont, or base resistance 

the cyclic shear strain along the pile to a critical shear strain for interface | 
_ skin friction y,,, and y,, would be expected to be on the order of 0.1 to 0.25 ; 
times y,. To estimate the cyclic shear strain y, in the soil adjacent to the 
> pile, use may be made of the analysis developed by Randolph and Wroth (15), 
which may be derived, for a pile in a deep elastic soil layer 
stad bed caphety of the pile, and will be 
in which p, = cyclic displacement of pile at a point on the pile shaft;d= ff 


pile diameter; v, = Poisson’ s ratio of soil; Le ‘embedded + pile length; vegme t" = 
In [sul — v,) L/d)]; and y= 1.0 if the soil is homogeneous and infinitely 
deep, or 0.5 if the soil has a a modulus which increases linearly with depth from 
' ‘zero at the surface. Eq. 4 will be only approximate if the soil ceases to be Ay 
elastic or pile-soil slip occurs, and in the latter case, may tend to overestimate ve 
a). the ‘cyclic shi shear strain of the soil. For the pile base, a similar approximate — 
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y skin friction Yep = scritical shear strain for skin 


For base resistance Yer =static failure strain’ for 
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4 FIG. : 2. —Basis of Total Stress Degradation Analysis: (a) + sae oe 


a "subsurface circular area. Assuming », 


in which | Pr = ‘= cyclic: displacement of pile base; a and d, = = diameter of pile 
‘The degradation carves Pow used for the ‘solutions described herein. 
are shown in Fig. 2(a), and have been derived from the data of Lee and Focht — 
mn Degradation of Soil Modulus £,. ~At least two approaches to » obtaining soil 
modulus degradation factors can be contemplated: (1) By making use of plots — 
of cyclic shear strain contours against cyclic stress level and number of cycles, 7 


such as those | presented by Andersen (1) for Drammen Clay; and (2) by using — = 


4 
Nodal 
| 
| 


in = ‘degradation parameter ‘depending 
on cycling axial strain. Because of the simplicity of the second approach, it 
= been used herein. Fig. 2(b) shows a plot of the degradation parameter 
e t versus cyclic strain ratio, , derived from the data presented by Idriss, et al. 
-. (8) for San Francisco Bay Mud It is assumed that the reference strain y. 
is equal to the static shear strain at failure, +,- In addition, because the = 
modulus is influenced by a relatively larger volume of soil than is the pile-soil oe 
skin friction, only a fraction w of the cyclic strain at the interface should be 
used in evaluating the modulus degradation factor. From an examination of © = 
4 results of finite element pre: ee a pile in a radially nonhomogeneous so 
obtained by Balaam, et al. (2), a value of w between 0.5 and 0.25 appear 
be applicable. For the solutions herein, w has been taken to be 0.4.0 
et Effect of Rate of Loading. —The undrained behavior of soils is influenced 


strength varies approximately linearly with the logarithm of time to failure. 

‘It seems reasonable to assume that the soil modulus, ultimate skin friction, 

and ultimate base resistance vary with loading rate in a similar fashion to undrained 

Strength. to allow for loading rate effects, factors 


= reference loading rate (e.g., for static load ane A= 
loading rate; and F, = rate factor, which the limited available data : suggest 
on the order of 0. 1 (range 0.0510 0.3), 


a -y typical offshore pile in clay has been analyzed, as shown in Fig. 3. The 
is soil parameters shown therein are considered to be representative of a normally- | 
consolidated clay deposit. In Fig. 3, degradation of ultimate skin friction as 
shown in Fig. 2(a) y,, = 0.02; degradation of ultimate base resistance as shown 
g Fig. 2(a) y, = 0.08; and soil ‘modulus degradation as shown in Fig. 2(b) 
y., = 0.08. The effects of variations in some of the pile and soil parameters 
will be examined subsequently. Attention is concentrated on the following — 
b quantities: (1) The ultimate load capacity after cycling, P,.; and (2) the cyclic 
“stiffness” (cyclic load divided by cyclic of the pile after N cycles, 
eo give the results s greater g generality, Po will be normalized with respect 
to P.,, the ultimate static load capacity of the pile, and K. will be nenutsined 
with respect to K .,, the cyclic st stiffness of the pile for the first cycle of loading. i 
‘The two other important parameters to be considered are the magnitude of a 
the “cyclic load P. (normalized as P./P,,), and the number of cycles N. In 
-allcases it will be assumed that the N load cycles will be of uniform neiteanginase eis 


(i.e., +P, to 
be considered: P,./P,, versus P./P,,; 


us? 


_P. range). Fig. 4 4 summarizes the definitions. of the ane 
cyclic response variables. In particular, the following four 1 ‘relationships will 


K./K., versus P./P 


us> 


_ All solutions have been obtained from a computer program CYCPL6, using 
10 shaft elements and one base element to divide the pile. For simplicity, the 


ultimate skin friction and base resistance in compression and tension are assumed ~ _ 
to be equal although in reality this would not be the case, particularly for 
the base resistance. Also, the ambient static load P, is taken as zero, i.e. ” 


the cyclic load has a zero mean value. 
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FIG. 3. Offshore Pile Analyzed: Standard Parameters 
re Bffect of Cyclic Load Level and Number of Cycles.—The dimensionless ultimate 


a joad and cyclic stiffness are plotted in Fig. 5 as functions of the number q 


Pile modulus Ep =21000 MPa 


“cycles N, for various dimensionless cyclic load levels P./P,,. Relatively little 
4 reduction i in ultimate load capacity occurs for cyclic load levels less than +50% 
of the ultimate static value, and data quoted by Bea, et al. (3) support this 
- conclusion. There is relatively little effect on cyclic stiffness, although as noted 
also by Boulon, et al. (4), the stiffness decreases continuously with epg x 
‘tm of cycles, even at relatively low cyclic load levels. However, for P./P,, 
= +0.6, failure during c cycling (P,. = P.) occurs after about 5000 cycles. - 

_ This sudden drop in ultimate load capacity with increasing cyclic load level é 
is further illustrated in Fig. 6, in which P,./P,, and K./K,., are plotted against z 
% /P., for various values of N. It is clearly seen that onde load | levels in 
excess of about 55% of P,, are likely to cause failure during cycling. _ ” a 

_ For N = 1000 ‘cycles, Fig. 7 plots distributions | of the degradation ‘factors, ~ 
cyclic load, - cyclic displacement, and cyclic shear stress , along the pile, for three 
values of P./P,,. As the cyclic load level increases, the degradation of both 


and E, becomes more marked, particularly near the top of the pile. Conse- 


4 


| ny, more load is then carried i ied in the the lower part of the pile. These characteristics _ 


4 
4 
| 
q 
| 
— 
g 

| 
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are also observed by Boulon, et al. (4) in their finite element analysis. Fig. 


clearly shows dt progress of “unzipping” of tl the pile (i. e., the loss of ultimate = 


(b) Pile ‘Settlement vs. Time 
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tc) Ultimate Load ‘Copacity vs. Time 
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4.— —Detinition Key cyclic Response \ Variables 
1 


\Values of 


» 


7 AIG. —§.—Effect of Cyclic Load Level and Number of Cycles on Cyclic Response 
dae a> the cyclic which ecoute the top. 
tics are found if the degradation factors and shear stress distribution are plotted — 
for i increasing numbers of cycles ata constant cyclic load level. Ole 
Effect of Critical Cycle Strain for Skin Friction, +,,. .—Fig. 8 ‘shows the effect 


j 
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ARY 1981 
the critical strain for skin friction on the ultimate pile capacity 
ond cyclic stiffness after 1,000 cycles. In all all cases, the critical strain for ultimate — x 
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- 6.—Ettect of Cyclic Load Level and Number of Cycles on Cyclic Response — 
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FIG. 1—tafluence of Load Level on Pile Behavior 
base resistance and soil modulus y, has been kept constant at a typical value | 
_ of 0.08. A decrease in y,, results in a significant decrease in both ‘ole 
us re and cyclic stiffness, with failure during cycling occurring at emailer 
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CYCLIC A) AXIAL \L RESPONSE 


load levels. The sensitivity thie i results to indicates that it would id 
\s be worth while performing mo model and field tests to obtain further data on values — 
4 of cyclic displacement necessary to cause significant degradation of skin friction, _ 
whence the value of y,,may be determined. = 
Effect of Distribution of Soil Modulus and Ultimate Skin Friction.—Calculations 
have been performed for both the standard case, in which E, and + .vary 
4 linearly with depth, and for the case of constant E, and r, along the pile shaft. 
The Static values of these parameters have been ‘taken as those at middepth 
of the pile in the ‘standard’ case, i.e., _ E, = = 50 MPa, +, = 0.125 MPa. Fig = 
_ 9 shows that the reduction of both ultimate load capacity and cyclic wanees ; 
| é with increasing numbers of cycles is much more severe for the case with constant — # 


| 


N=1000 


02 04 06 08 10 


8. —Etect of Critical Shear Strain for Skin Frtetion, Ves 


—Effect of Distribution of Soil Modulus. and | Skin F Friction 
_ £E, and +,, and in fact, failure in this case will occur after about 16000 cycles ‘ 
it P ./ P,, = +0.5. The implications of Fig. 9 are of some practical importance, __ 
as it is ‘clear that piles deriving significant frictional support from near-surface a by 
_ soil layers under static loading conditions may lose a significant amount of a 7 
this support due to the severe cyclic degradation which occurs near the top 
‘ of the pile. Piles in overconsolidated clays may fall into this category, a 
Effect of Loading Rate.—The preceding solutions assume no effects of the 
ae rate on ultimate skin friction or soil ‘modulus. An indication of the 


4 

4 
>. 


‘jnfiuence ¢ rate” effects is given in Fig . 10, in which quantities, in 

dimensionless form, are plotted against a dimensionless loading rate factor, — 

_ P.,/T,, in which T = cyclic loading period, and \, = reference rate of loading 

for which the values of static soil modulus and skin friction are determined. il 

_ As would be expected, an increase in the rate of cyclic load application (and = 
hence in P. s/Th, ) significantly increases the vu ultimate pile capacity, although 
-y does not have much effect on the cyclic stiffness ratio K.,/ K.,. The beneficial 

effect of increased loading rate becomes more pronounced as P | Pus increases. 

The increase in pile capacity with increasing load xate is consistent with data 


and Matlock (7) have found considerable differences between the p mer a ol 
of piles subjected to ‘‘one-way’’ and ‘“‘two-way’’ cyclic loading, with cyclic 
_ degradation being - much more severe in the latter case. One of the reasons _ 

for this behavior i is that larger cyclic stresses and strains are developed i in two- put 


Sol'ns. for no rate effects 
— 
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LR), = 


(LR), = Reference rate of leading 


& FIG. 10.—Effect of Loading | Rate 


skin friction and modulus will occur. Thus, the same maximum load 


is dependent on 1 cyclic load level, the ultimate load capacity and cyclic stiffness 
_ for one-way loading for a maximum load P,,,, will in fact be the same as those - 


max 


for two-way loading with a maximum load of 0.5 P,,,,. 


A series of small-scale model | pile tests was carried out in four different — 
soils: Sydney Kaolin (LL = 47, PL = 18); Hurstville Clay (LL = 65, PL 
> _ = 20); Redfern Clay (LL = 42, PL = 22); and Darlington Clay (LL = 55, 
PL = 33). In each case, the soil was remolded and placed in a cylindrical - 
152 mm (6 diam mm in.) deep. The piles were 


& 
i 
g 
| 
4 
dee’ 
may ve reac in DOLD Cases, 1055 OF UllimMate load Capacily and pile sti — 
| 


‘so 20 9 mm (0 78 in.) diam and 250 mm (10 in.) long, and were = 
jacked through the clay so that the end of the pile protruded from the bottom __ 
Of the soil. Consequently, only the skin friction characteristics were measured. _ 
Fe After allowing a suitable rest period (typically two to three hours), the following __ 
procedure was carried out: 
1. The pile was loaded to failure in compression over a period o of about — 
2 min in order to obtain the Static ultimate load capacity P., (failure being 
defined as a deflection of about 2.5% of the diameter), 
* After unloading and leaving the pile to rest for | hr, a cyclic load a 
pan 


Hurstville ‘Clay 


Kaolin 

= 


_ Theoretical “Curves (see Table 1 for 


after 
half was applied to. the by a piston operated 
b a compressed air-actuated relay valve; 1000 cycles of load were applied, the 
average period of loading being2.5sece | 
4 3. The pile was loaded to failure in compression, generally over a period 
. of 2 min, in order to obtain the postcycling ultimate load capacity P,... In the” 
case of the pile in Darlington Clay, the time to failure was typically about © 


_ Because of the nature of the apparatus, the mean load about which cycling — 

occurred was about 37 N (8.34 Ib) (approximately 10% of the ultimate static 
load). The load was measured using a strain-gaged proving ring, and the — 
\ deflection of the pile head by a displacement transducer. Load versus time > 
- deflection versus time relationships were plotted directly ona chart recorder, 


| 
| 
‘4 


while load- plotted on an x. Y 
Fig. 11 summarizes the data on ultimate load capacity after, 
cycling, and plots P,./P,, versus P./P,, for N = 1000 cycles. The experimental 
‘Tesults for the four ‘soils are reasonably consistent and indicate little reduction A ; 
in ultimate load capacity unless the half amplitude of cyclic load P. exceeds _ 
about 60%-70% of the ultimate static load capacity P,,. The results for Darlington 
Clay, indicating values of P,./P,, greater than unity, may reflect the influence — 
of the higher loading rate used in these tests. § 
‘Fig . 12 shows experimental results for the variation of cyclic stiffness with .* 
- number of cycles, for a cyclic load level of f approximately +0. 45 Pus: agp = 


there is 


ar vf he reasons 
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Theoretical Curves for 70.42 


12.—Experimental and ‘Theoretical Results for Cyclic Stiffness Versus Number 


of Cycles 
_ The influence of cyclic load level on cyclic pile stiffness is plotted in Fig. — 
13 for N = 1,000. As would be expected, the decrease in cyclic stiffness (as 


‘ expressed by the ratio K./K’) becomes more pronounced as the ‘cyclic load 

Because of the difficulty of assigning accurate values to the various parameters | 
required in the analysis, theoretical solutions have been obtained for a range — 

of parameters. The soil modulus E, and the static ultimate skin friction 1, i 
e assumed constant with depth, the values being chosen to give r reasonable 

used are shown in Table 1. The degradation of skin friction and soil modulus _ 

has been determined from Fig.-2, 
‘The theoretical solutions thus obtained are shown in Figs. 11 through 13 
S show similar trends to the experimental data. Curve 5, which includes © 


allowance loading: rate effects, is relevant only to the tests on Darlington 
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1,000 cycles, 4 
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The lowing points are worthy of note 


1. The theory predicts a very sudden in ultimate load capacity P uc after 


TABLE 1. 1 .—Parameters Used in in in Theoretical Solutions for Model ~ as . 
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FIG. 13.—Effect of Cyclic Load Level on Cyetie P Pile ‘Stiffness; N = 1,000 Cycles HO = 
2 The value of y,, (the critical shear strain for strain friction) to obtain 


Jt 


agreement with 1 the experimental results is 5 quite : small, on the order r of ( 0.01. 0.01. ae 
ae the “uncertainties still existing in the values ; of the input parameters, 


| 
: 
Proportion of 
| Curve 
; 
a 
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a egfoement that in the ‘model tests. It is interesting to note 
the present total stress solutions are in better agreement with the test results _ 
than are the effective stress solutions presented previously by the writer (13) 
when pct ge the results of the tests in Darlington Clay puede 


‘The analysis presented herein is an extension of an elastic- ‘based analysis" 


effects of cyclic loading, and which involves an iterative analysis to determine : 
the ultimate load capacity and the cyclic” stiffness of the pile. after | a given 
number of load cycles of a given ‘magnitude. 
_ Some solutions have been obtained for a typical offshore pile, and f and from 
the following observations have been made: 
— ultimate load capacity and cyclic stiffness decrease with increasing 
4 “number of cycles N and increasing cyclic load level. The effect of N becomes - - 
; E more significant when the cyclic load P. (half peak to peak) approaches 50% : 


.. 2. Cyclic degradation begins at the top of the pile and progresses downward 


as P. or N increases, ‘Tesulting in gradual transfer of load to lower parts” of ai 
3. The crucial factor in determining ‘the amount of cyclic degradation i is the 
critical shear strain for skin friction, y,,. ‘there is little data 


— the many assumptions and the simplified nature of the analysis, . 
- does give some insight into the behavior of cyclically loaded piles, and = 
behavior which consistent with the results of small-scale model tests 


further data | become available on some of the key input ‘parameters regeired, 
“~~ may provide a basis for assessing cyclic load effects on off fshore piles. — 
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i 
' D, = degradation factor for ultimate base resistance; mae 


K, = cyclic stiffness of pile; 
Ka = cyclic stiffness of pile after | cycle; 


P. ultimate static load oad capacity; 
mex = maximumiosd; 


= = value of p,, after cycling; » 


cyclic loading period; 


proportion of interface cyclic stress or me used for modulus 


critical value of shear strain; bat 
~y ‘reference rate of loading for which : static parameters are determin 
te rat static shear strain for failure of triaxial sample; : 
= Cfitical shear strain for skin friction; 


% = value of +, after cycling; 


= static falue of and < A 
y= number defined in in Eq. 4. 


a 4 


| § = degradation factor for soil modulus; 
E. = soil modulus after cyclicloading; = | 
| 
4 

| 
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EMBEDDED ANCHOR RESPONSE TO LOADING 


‘ Andreas Andreadis, ' Roger C. Harvey, and Eidos 

ij 

Embedded anchors or ground achors whichm may consist of a deeply embedded 

- plate, bulb, or flukes connected to the anchored structure by means of a cable 


with both research and practice (8,10,19 21). ‘Their potential as an acceptable _ 
anchor for close station keeping in the marine environment has also been ; 
recognized and a number of systems have been | developed (6,14,17,26) which, 


as ocean operations move into yet deeper waters, may prove increasing importance 
' _ As a sea bed anchor can be subjected to a wide range of loading a a 
including sustained or static, repeated, and sustained-repeated combinations, 
A it is essential to obtain a thorough understanding of soil-anchor interaction. 
‘In particular, an appreciation of the life of the system and influence of installation 
(3) are important factors as the ability to efficiently design these 
anchors may overall economics of future facilities. 


The anchor “one on in the investigation are shown in ‘Figs. 
Those in Fig. 1 were designed to be installed vertically using a fluidization | 


process. Some cylindrical anchors contained pore pressure gages and a system _ 


_ capable, if required, of eliminating cyclic creep effects > ee ar 
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Fig. 2 shows a model was or (with closed flukes) 
ay driving. At a predetermined depth the flukes were unlocked by activating _ 
an explosive bolt. The final phase of driving caused the soil bearing pressure a 

it to force the flukes open until the | contained angle was 90°; the locked flukes : 
then offered a substantial withdrawal resistance during subsequent testing. 

_ The tests were performed in a steel tank measuring 8 ft (2.44 m) in diam > 
and 4 ft (1.22 m) deep, positioned on a vibrating table and filled with a fine- medium, 

uniform, saturated sand (Borough Green sand) with engineering 


Re varie in Table 1. Additional tests were carried on in Bush Farm 


ee 


Merc} fue 
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— F FIG. 1 —Cylindrical, Conical, and Plate Anchors 


also shown in n Table The table been drawn up in “accordance with the 
‘Unified Soil Classification, with the visual information following the Rittenhouse 

and Krumbein methods. The standard methods used for determining the soil 


properties were in accordance with the American Society for Testing and Materials _ 


The size of the test container was chosen to ) reduce boundary effects 
: poe been shown (22,27) (for deep anchors embedded | in dense sands) to be 
present at a distance in excess of 10 anchor diameters from the loaded anchor. 
ee method of soil bed preparation had to take into account the 
4 large masses of soil and water involved, a total of approximately 10 metric 
: Repeatable soil beds of known properties were obtained using vibeshocy, 
2, 5) after installation: of the model an chors. . Soil pressure gages 
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a = = EMBEDDED ANCHOR RESPONSE 
show thal tbe ‘ners with 

q The complete testing equipment is shown in Fig. 3 and consists of the test 


tank mounted on a vibrating table with a servo- -controlled hydraulic loading 


i = These tests were carried out to determine the influence of anchor eles 

other factors involved on the ‘Static Capacity and response 


1g 2.—Fluke | Anchor and Explosive Bolt 


to different levels of static loading. A pull-out Tate of 0.5 mm/min, “Te 


_ Influence of Relative Depth and Anchor Size.— —Fig. 4 presents typical static 
test data in terms of the static breakout factor, N,, against the relative depth, ie 
A, (= in the anchor and = the lateral dimension 


atlacne the rigid steel Irame. applied, by means of a Tunclion 
, 
= 
| 

| 
| 


ey, 


the ‘ultimate “uplift ‘capacity increases sat a rate greater then in proportion to 

D’, but as the depth of embedment i increases and failure becomes more localized _ ~e ~ 
there is not the same dependence on depth as for the shallow anchor, and 

: the ultimate load-relative depth relationships tend to become linear, 2 3 
. ce has been established (27) that significant scale effects are present when 

sd comparing small-scale anchor results with medium and larger scale results. “a 

Although breakout factors increased ‘significantly for anchor plate ‘diameters 

Tess than 1.97 in. (S50 mm), it appeared that scale effects were considerably fx 

reduced, with breakout factors substantially unchanged, for diameters greater 

than 1.97 in. (50 mm). This was endorsed by field tests carried out on prototype 4 

anchor flukes embedded in dense sand in which good correlation was obtained  __ 

4 between model and prototype b breakout factors (2) (Scale factor = 3. 3. 4). Similar i 

results obtained by other authors (7,18,25 28) who suggest that cautious compari- ‘ 

a son between medium- scale model laboratory tests and larger scale field tests” . 


Maximum dry density 1.80 g/cm’ 
Minimum dry 1.44 g/cm’ 

Permeability (at D, = 10%) 2x 10°? cm/s 3.5 x cm/s 


Fig. 5 compares the test results with a theoretical predictions 
of the ultimate static uplift capacity. Uplift prediction for shallow anchors ee 
i not, in general, yield good results mainly due to the assumption of a unique 
“3 failure surface shape for all soils, in addition to an inadequate representation — 

2 the soil forces at failure. The complex soil-anchor interaction problem ; appears $ 
to be oversimplified by these methods. Meyerhof and Adams’ (22) deep anchor _ 


_ analysis yields comparatively good results although at greater bmg the — 


4 
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‘Maddocks: (20) that the behavior | of a plate anchor in 


same immediately above the with the lifting of large zones of sand above 


scements at failure i ‘is miso 


| 
; Jecic’c (? a H ma 
J 
3 predicted by the expanding cavity method (11,29,30). Kupferman’s (18) solution a 
takes into account the observed failure mechan 
i. - much larger zone of soil observed to suffer displ 
. Examination of the influence of relative dept 
- 8 the movement due to static loading and movement A 


RRESPONSE 
at static failure, respectively, show that the relative movements increase with — 


_ relative depth in a manner similar to the increase of the breakout factors. A 


g Equipment 
as 


; hor load-displacement characteristics, shown in 
_ Fig. 6, suggests that within the limits of experimental accuracy, a unique — 
relationship for cylindrical anchors embedded in undisturbed dense sands exists 


| 
| 
| 
i 
2, 
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of the soil and schohintatly on the embedment method in addition to the total 
a cyclic movement the anchor has experienced, as will be later explained. : 
Influence of Anchor Shape. -—Static | loading tests were performed w with a number 
a anchors and typical results are shown in Fig. 7. The increase in 1 uplift « capacity 
a increase in plate thickness can be attributed to a reduced movement of # 


Cyingricel Archers 

B= 75mm Or = 66% 
Or= 66% 


eB: 


Conical Model 
* B=75mm Or = 66% 
Plate anchor B= Dr = 66% 
Fluke anchor 8 =140mm Dr= 66% 


the early stages of loading the: anchor exhibits a . similar load- -displacement ehevior 
to the plate and cylindrical anchor. However, as loading increases, a sudden 


ultimate capacity. This reduction is considered to be due to the increased flow 
of soil beneath the anchor, thus. reducing the effective confining soil pressures. i d 
_ Comparison of the cylindrical and fluke anchor breakout factors in Fig. 4 : 
shows that the fluke anchor displays relatively higher breakout factors at shallower 
oF depths. However, as the relative depth increases, the more ‘‘aerodynamic’ g 
wedge | fluke shape apparently allows local shear failure — to more easily ‘ 


_ 
ion in the Ra city he noch nes = 


develop at reduced This ‘conclusion w was endorsed comparison 
& the post-peak behavior of the anchors at shallow depths in which the fluke — 
anchor displays a deeper failure mechanism than the cylindrical anchor. Clearly, — 4 
4 as the apex angle of the fluke wedge reduces, the transition point between — ae 
shallow and deep failure will move closer to the soil surface ee Bea 7 
Be _ Extent of Soil Disturbance during Testing. .—Electronic soil pressure gages were 4 
embedded at specified locations in the soil mass ss in order to monitor the soil ¥ 
Jp : 


Cylindrical model anchor 
Cylindrical model anchor 
bag Theories 
Mors (1959) 
Balla (1961) | 
Mariupol’skii 
 Meyerhof & Adams (1968) 
Vesic (1971) 


_evione Relative depth A= -0/8 


FIG. 5.—Comparison of Static Load with Wok 
+ Fig. 8 shows the distribution of horizontal stresses at the tank wall during — 
_ the testing of a 4.92 in. (125 mm) diam cylindrical model anchor. This is equivalent 
to 9.6 times the anchor body diameter from the tested anchor and similar results © 
were obtained during the testing of a 2.95 in. (75 mm) diam anchor. It is evident 
that at such distances significant horizontal stresses are generated in the soil 
mass. The stresses increase rapidly from zero at the soil surface to a maximum eS 
~ value and then decrease again to zero just below the level of embedment of a 
the anchor plate. The position of the peak value depends on the distance of 


a o monitored stress profile from the loaded anchor, e.g., at distances greater \ 
than five anchor diameters, ‘maximum horizontal stresses a depth 


= 


aa to half the embedment anit of the anchor. Closer to > the loaded nce 
Kd body horizontal soil stresses substantially increase (see Fig. 9) and the position — 
_ of the peak value shifts nearer to the anchor body. The pear nt of the — 
= pab brs 
age 


Q/Quit 


Al depts 
Both tested sands Dr = 


a7 
movement 


Bt = Anchor movement to failure Re 


FIG. —Stati Load 
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FIG. —Static Load Tests: Uplift Load- Anchor Movement Relationship 
- appears to be significantly affected by both anchor size and Telative depth o 


‘Fig. 10 shows typical increases in the horizontal stress at a “specified location 


iy 
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_ EMBEDDED ANCHOR RESPONSE 
in the soil mass a static loading test. Horizontal soil stresses measured 
as distant as 16 anchor diameters from the loaded anchor continue to increase _ i: 
in a mode characteristically similar to the load-displacement behavior of the 
anchor. After m= anchor resistance — and en takes place at constant: 
B = 125mm “hee 


& 


Wit 
Q/Quit 
5S 1 15 2 25 30 


B=125mm —o— 
B= 75mm A=8 
Borough green sand 

Relative density Or 
Stresses measured at A=4 
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FIG. 9.—Static Load Tests: Horizontal Stress- Relative Distance fom Anchor oe Ra; “s 


res 
“load, the soil stresses continue to increase at a reduced rate and tend to : a 


constant value when the anchor resistance falls. dike 
4 This behavior suggests that a much vom mass of sand is affected by mens, ? 


8.—Static oad lests: Horizontal Stress Profile on Tank Wall 
a 


_ anchor plate (7,20) are probably unable to reveal the entire zone of stressed 
soil, as dense sands appear capable of mobilizing and transferring substantial — 
stresses without easily detectable grain movements. For anchors embedded at 
: relative depths, less than eight horizontal stresses become rather insignificant _ 
_ at distances greater than 15 anchor diameters (see Fig. 9). Deeper anchors appear 
to develop significant horizontal stresses over a much larger soil zone and 
10 


ess 


Horizontal str 


10.—Static Load Tests: Horizontal Stress-Relative Anchor Movement Relation 
indicat model anchor 
— unloading 


re water 


“we 


Anchor Loading-Unloadin 


this zone is restrained some interference 
anchor could be anticipated (27), 
__ Vertical soil stresses, developed due to anchor testing, ' 
These stresses were, as expected, substantial above and near the anchor body. 
_ However, their variation with distance from the anchor showed a very marked | 7 
reduction compared to the horizontal stresses, and at distances greater than — 


five anchor diameters no significant effects could be 


* 
used to obtain accurate representations of sand grain displacements around the =] 
he 
| 
| 
Fy 
FIG. 
a 
with 
FIG. 11.—Pore Water Pressure during Fast 


_ EMBEDDED ANCHOR RESPONSE 
Influence of Pore Water Pressures.—The selected testing rates were ae 


to provide essentially drained tests. . However, if stresses are rapidly applied 
to o dense cohesionless soils, only a limited dissipation of induced negative pore — 
pressures will take place resulting in partially drained conditions, and therefore 
increased initial holding capacities; 
_ Fig. 11 presents typical readings obtained by the pore pressure gages fitted — 4 
below the anchor plate edge during a very fast load application (snap load). j 
_ Negative pore water pressures appear to o accumulate rapidly, reach a maximum — 
- Maximum snap loads attained by the tested anchors were more than double. 
their static drained capacities. Sudden relief of the load resulted in a rebound 
of the overlying soil mass associated with a decrease in the volume of i oe 
voids and a generation of positive pore water pressures. Induced negative or 
ositive pore water pressures appear to dissipate at at exactly the same rate (Fig. 
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Diameter B= 75mm 
Relative depth A-= = 

Borough green sand ua 
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yclic movement 4 
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S FIG. 1 2.—Repeated Load Tests: Relative Anchor Movement (Log Number of Cycles | 

_ As medium dense cohesioniess soils dilate markedly less than denser soils 

_ whens subjected to sh shear s stresses, increases in undrained anchor holding capacities 

 peitiieel oe pore pressures, and the undrained holding capacities i are e expected to 


be reduced in comparison with the drained case. 


_ Each anchor was subjected to repeated loading which consisted ofa contiontes = 
- sinusoidal loading pattern. Various dif! ferent einen loading and sustained- 


| 
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: 
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JANUARY 
repeated loading combinations were used and similar results were found for 
_ standard cyclic period of 2 sec was adopted as a result of preliminary 
investigations employing several loading periods from 1-30 sec. It was judged 
_ that a frequency of 0.5 Hz introduced negligible inertia forces and testing was ‘ 
therefore essentially nondynamic. Additionally, the pore pressure gage readings _ 
aaa drained conditions. The adoption of a 2-sec loading period advanta- 
geously permitted accelerated testing and approximately 500,000 loading cycles : 


rf the repeated loading « characteristics of shallow an ould deeply « embedded fluke, | plate, 
and cylindrical anchors (1,18,20,23,24) suggest that the basic response of each 
= 60% 50% 40% Cylindrical model anchor 
L; Diameter = 75mm 

Relative depth ae | 

Borough green sand 

Relative density D, = 66% 
Repeated loading 0 = 0. 


cyclic movement, 


= an 
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0° 


FIG. Load Tests: Relative Anchor Movement == Rate of Relative 


anchor is similar. In the early ‘stages: ‘feet the upward 
¢ accumulates rapidly. However, if cyclic loading amplitudes are kept below a 
critical value, the rate of displacement of the anchor continues to reduce ell 


never ceases. For shallow anchors, D/ B < 8, an increase in the rate of movement | 
‘may eventually occur (18) due to reduction of the relative depth as the anchor : q ' 
displaces upwards. For very deep anchors, D/B > 40, however, the soil surface — 7 
_ influence is limited and as the test progresses the system appears to become __ 
' more stable (9). On the other hand, if cyclic loading amplitudes exceed a critical — 
value, the net upward displacement for successive load increases, 


- Fig. 12 presents a typical repeated loading response of a cylindrical model 


anchor embedded in dense Borough Green sand. The data is presented in terms 
of the accumulating relative anchor movement, Ad, against the logarithm of 
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"EMBEDDED A ANCHOR RESPONSE 
- number o of load cy cycles, log N. The an anchor is ; subjected ‘to loadings which an are 
in the form of sine waves op. between zero load and specified peak | . 
relative cyclic load values, Q./Q,,,, expressed as a percentage of the ultimate _ 
drained static capacity of the anchor. In all cases, an unceasing anchor movement _ 

develops, its rate initially tending to decrease with time. As the test progresses, 
-_ however, the gradients of the curves increase indicating approaching failure. < 
- Repeated loading behavior of a particular anchor system can be represented ; 
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in a, ‘and c= cyclic constants. 

Cylindrical model anchor its 

Diameter B= 75mm 

Relative depth $i A= 8 oF 

Borough nd 

Relative density Dr = 
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FIG. 14.—Life of Anchor for Different Relative h Movement Lit Limits as Function « of vf Guaie 
Load Amplitude = 


a Similar families of f hyperbolic curves were found to exist for different relative — 
* depths, sand densities, anchor sizes and shapes, installation methods, sustained- > 
repeated loading combinations, as well as for varying degrees of soil flow around — 
and under the anchor. It appears, therefore, that the proposed method of analysis | 
: could be formalized to express the basic cyclic loading displacement and failure — 
characteristics of the majority of embedded sea bed anchors. © regia! aomoosd 
; e Fig. 13 shows the. accumulated relative anchor movement plotted against the 
system, the definition of which has caused difficulties in the past, can be defined — 
_as the turning point on the rate of relative movement against the relative anchor _ 
‘displacement curve. Past the critical point, the net anchor movements per cycle a 
_ increase and failure takes place rapidly. 


: strain to the minimum point is the “critical cyclic rela relative 


4s of relative deformation on a semilogarithmic scale. Failure of the anchor 


— 


“movement, and a appears tov vary both with 1 depth of embedment relative 
p cyclic load. For design purposes, however, a mean value of Ad. equal to 0.1 
: -. the tested anchors can be used. Accumulated relative anchor movements | 
excess of the critical value will imply imminent failure. Subsequent testing 
= that anchors experiencing total strains below the critical value display — 
negligible drops in their ultimate static capacities. dt 
The definition of repeated loading failure clearly implies that cyclic strength ee 
deterioration can be best explained as a strain dependent phenomenon. Further 
proof was obtained by os post-cyclic response to loading, in addition to subjecting oe 
the anchors to groups of cycles with the same frequency but with varying A 
stress s amplitudes. It appears that the deterioration i in n the es loading- relative . 
= 


‘movement, whether developed by a few ‘strong or - many small stress “pulses. 


“Fig, 14 shows the progressive accumulative movement of the anchor as it dy. 
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FIG. 15. —Repeated Load Tests: Hysteresis Loops 


"developed after three simulated consecutive two- y-day storms with equivalent 
- uniform repeated load amplitudes, Q./Q,,,, of 20%, 40%, and 20%, respectively. 
During testing it was observed that the area of the hysteresis loop decreases 


in size (see Fig. 15). At near stable shape of hysteresis loop i is } achieved minnll 


continues to accumulate. After the critical relative 1 movement (Ar di is reached, Bs 
_ the permanent upward anchor displacements increase and the hysteresis loop 
becomes larger until the anchor is pulled from the soil, 
a The repeated loading characteristics of an embedded anchor are significantly _ a 
i affe *d by the installation depth and the upward movements of an anchor _ 
a — subje ted to a given relative cyclic load increase as the relative depth i increases. P 
Although it appears that a weaker repeated loading behavior is displayed by 
_ deeper anchors, a different picture of anchor efficiency emerges if repeated <_< 
loading characteristics are compared in terms of absolute rather than relative 4 
eye loading amplitudes. Thus d changes between zero and OkN correspond 


= 
q 
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q GT1 
ma to | 60% of the: ‘uplift capacity of an anchor ‘for which B = 2: 95 i in. (75 mm) 
and d= 8, but only 23% of the uplift capacity of the same anchor installed 
ee at \ = 12. The former, when cycled at Q./Q,,, = 60%, will fail after 35,000 
load repetitions whereas the latter when cycled at Q. / Qu = 23% will fail at 
175,000 cycles, a a five-fold increase in life. Greater improvements in life. ne 
- favor of deeper systems were e observed when the rep repeated loading characteristics 
a of anchors embedded at relative depths of 4and 8 were examined. 
_ Forthe tested anchors the critical cyclic relative movement, Ad., imaepacindl 
to failure could be regarded as 0.1. Tests show however, ‘that values 
‘diminishes, large cyclic relative displacements 1 be attained before the rate 
es of anchor movement starts to accelerate. This conclusion has been endorsed 
_ by other writers (9) who have shown that no increase of displacements per 
_ cycle were monitored for very deep anchors (r = 87) although anchor plate _ 
movements were inexcessof0.5B. = = | 
bes Although the basic mode of response is essentially unchanged, | different size 
repeated loading behavior. ‘The gbeerved scale ‘effects can be attributed, in 
part, to the cyclic creep phenomenon, i.e., the flow of soil around the anchor 
Z a _ tending to partially fill the cavity formed under the anchor as the test progresses. — 
_ ‘The rate of soil flow is largely dictated by the relative elastic displacement — 
' per loading cycle which is greater for the larger anchor system. _ Therefore, 
= a faster cyclic strength deterioration | for this anchor must be expected. _ Experi- | 


mental observations with an anchor ‘system eliminating cyclic creep (4) verifies 
_ Scale effects were also believed to. arise due to differences in the peak cyclic - 
_ pressure developed on the anchor plate and to the increased grain ag ve el 
“taking place above the larger anchor body. Subsequent testing with an anchor 
“model eliminating cyclic. creep suggested that, provided the soil flow is eliminated 
the cyclic load amplitude is opt low, anchor behavior is essentially 
unaffected by size, 
Anchors subjected to sustained- -repeated nation are subjected to = bel 
«Mascon» at the commencement of the test when compared to those subjected _ 
to only the repeated loading component. This weaker behavior is considered — 
to be caused in part by the larger cavity created under the anchor body and 
s by the higher imposed peak loading. As the test progresses, the stress _ 
: system developed in the sand above the anchor body is such as to reduce 
: _ the rate of displacement to levels experienced by the same anchor system subjected oe 
Bes to the repeated loading component. Thus the cyclic relative movement ‘ 
Of | anchors subjected to sustained-repeated loading combinations, AA, | can be be 
derived » reasonable from the ¢ equation oa 
in which Ad = the cyclic relative movement of the anchor subjected only to 


the repeated loading component; and d= corresponding cyclic loading 
It is interesting that the relative elastic displacement per loading cycle a a 
asthe sustained load amplitude increases and the most severe elastic displacements 
rf occur when the anchor is completely unloaded ties each cycle of loading. 


| 


i solely to repeated loading could display a marginally weaker response. Similar __ 
have been reported by other workers (20,9), 
Influence of Anchor Shape.—Comparative tests were performed on cylindrical * 
and conical anchors and identical repeated loading behavior was found. In the ; 
J “case of the fluke anchor, response was similar except that the fluke anchor 
It is of interest to record that Kuptinnes ( (18) performed repeated loading 
bf - on model fluke anchors embedded in dense saturated sands. A direct — 
comparison with the information obtained during the present work is 7 
in Fig. 16. Considering the differences in fluke size and geometry together — 
with the variations in the soil properties, preparation, and testing procedures, — 
the overall correspondence appears to be reasonable. 
Fluke anchor models 
B=140mm,A=6 on 
B=76mm, A=4 (after Kupferman, 197) 


45° Or = =88% 


B= 76mm, A=4 lafter Kupferman, 1974 ) 
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. 16.—Repeated Load Tests: Comparison with Kupferman’s Data (1974) 
Stress Distribution around 
- pattern identical in form to, and of the s same e frequency a: as, , the loading Gans. 
Their magnitude is mainly related to the amplitude of the loading change during» 
each cycle and it increases substantially as the loaded anchor plate is approached. 
a A comparison of the static and cyclic horizontal stresses that develop in 
the soil mass when the anchor i is loaded from 0-40% of its ultimate uplift capacity — 
“either once or in a repetitive manner is shown in Fig. 17. It is clear that after 
‘the initial application of load the effects of the repeated applications of loading — 
are confined to a limited soil zone surrounding the anchor body. 
_ Maddocks (20) carried out stereo-photogrammetric observeiions of plane strain - 
anchor models and concluded that sand deformations associated with stresses 


developed by each cycle of loading were me small and restricted to 


a 
10 102 10? 10° 10° 10 
— a limited 


— Finally, the vertical cyclic stresses dovdaged ts in the s soil during repeated loadi * P 
showed a very marked reduction when compared to the corresponding vertical _ 
Post Cyclic Response to Loading.—Tests were pr performed i in which the | anchor, 
- after having been subjected to repeated loading, was then loaded at a constant 
rate to static failure. A stiffening in the anchor displacement characteristics | 
was found, as compared to that demonstrated by only static loading, which — 
was associated with the soil density changes and supported by evidence of © 
a me reductions in the area of the hysteresis loops. Anchors subjected 
* to a few strong or many small stress pulses, or to any sustained-repeated load | 
, combination, increase their resistance to displacement at a decreasing rate as 
the cyclic movement accumulates. They attain an ultimate hardening value when _ 
_@ particular cyclic strain limit is achieved. Past this limit no further stiffening 
Deol 


ini appears that this cyclic strain limit can be related to the egervet 4 
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eee horizontal stresses 
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17.—Comparison of Static: Horizon tal Stre 


- relative movement as corresponding to anchor failure. Indeed, the inability of 
the ‘sand mass after failure to competently resist further applications of load, 
-_ denoted by gradual increase in the rate of cyclic movement, appears to 3 
4 associated with its inability to mobilize greater stiffness. Anchors deeply embed- 
ded displace more before failure occurs, thus increases in their stiffness 


4 
_, 


lic 
| 
— 
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‘This work has attempted to provide a better understanding of the behavior _ : 
_ ¢ of embedded anchors. The anchors were subjected to both sustained and 
Sustained-repeated loading combinations. Attention has been focy 
of important features relating to both general respons 
following brief conclusions may be considered applica 
| 


The mode de of fa failure of an anchor subjected to static loading is mainly es. 
le by the relative depth of embedment, soil relative density, and anchor 
Static breakout factors and relative anchor movements to failure increase 
_ with relative depth of embedment at shallow depths, tending to an 
approximate constant at greater depths, ga 
_ Significant horizontal stresses pa apg in ‘the sand mass in a zone which — 
can extend considerably in in excess 0 al 10 anchor diameters either side of a ‘a deeply 
2. The repeated applications of load to an anchor substantially reduce its 
resistance to load and a nonrecoverable movement develops. A realistic failure 
criterion can be interpreted using a rate of deformation technique which is 
defined by unacceptably high relative relative movements of the anchor resulting * 
in significant reductions in post-cyclic pull-out capacity. brite 
_ Anchor movement, , due to successive repeated applications of load at varying © 
7 amplitudes, can be predicted by using an experimentally observed strain dependent _ 
__ Stress measurements have shown that after the initial application of the load, 
the effects on the soil caused by the repeated applications of the load are — 


confined to a limited zone surrounding the body o of the anchor. + ee eee 


It has been important to to investigate response to both sustained and repeated _ 
loading, not only to gain a more complete appreciation of behavior, but also 
_ because in- eereyes gee under marine conditions may be sustained only, e.g., 


for the design of deep water moorings necessary for the future offshore industry. 
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‘There ar are ¢ several areas of California w where subsiding s soils are of signifi icance. a 
x Subsidence may occur because of withdrawal of fluids (usually deep subsidence), j “il 
_ the drainage of peat lands, or by the application of water to moisture- deficient, 
low-density soils (usually near-surface subsidence) (7). Only one type of — : 
dence, i.e., subsidence caused by. the compaction due to wetting, is analyzed 
in this paper. This process is often referred to as hydrocompaction and has 
widespread subsidence of land surfaces (5). 
__ Deposits that subside because of hydrocompaction are generally one of two 
7 types: (1) Loose, moisture-deficient alluvial deposits; or (2) moisture- deficient 4 
loess and related eolian deposits. These deposits, laid down with an open structure 
and high porosity and buried with a moisture deficiency ands: sufficient intergranular 
strength to support an increasing overburden load, are readily susceptible to 
_ compaction when wetting reduces their dry strength. They are generally found 
in regions where rainfall seldom penetrates below the root zone, and they 
. invariably occur above the water table. Unlike land subsidence resulting from 
compaction at depth- -of-water- -bearing deposits: due to a reduction of fluid 


~ hydrocompaction i is a near-surface phenomenon. Since w wetting progresses from 


_ progresses downward with the advancing waterfront. This near-surface subsi- 
_ dence is a serious problem and has damaged canals, roads, pipelines, and 
transmission t towers in act areas. It presents particularly serious problems 
i : nai 
in the ction aint ce of large canals. mea 
x 


TChf., Geotechnical Branch, Div. of Research, “United States Dept. of the Interior, 


Pi the land surface downward, hydrocompaction begins near the surface and — 


Building 67, Denver Federal Center, Denver, Colo. 80225. 


Note.—Discussion open until June 1, 1981. To extend ‘the closing date one 
3 written request must be filed with the Manager of Technical and Professional a 


5 Water and Power Resources Service, Engrg. and Research Center, P. (0. Box 25007, 


ASCE. Manuscript was submitted for review for possible publication on June 24, 1980. 
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3 The San Luis Canal is part of ” “California Aqueduct (see Fig. 1). The 
canal was designed and constructed by the Water and Power Resources Service _ 

- a joint venture with the state of California. (On November, 6, 1979, the ane 
Bureau of Reclamation was renamed the Water and atiiead Resources Services 
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FIG. 1.—Location of Prominent Features, Central Valley, California 
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‘The text of this publication was prepared lame; 
all references to the Bureau of Reclamatio re to 


be considered with the Water atl Resources Service. ) The 
San Luis Canal is 107 miles (172 km) long and has a capacity ‘Fanging from cae a 
100 cu ft/s (370 m’*/sec) at the headworks to 7,750 cu ft/s (220 m? m’*/sec) 
at the lower end. It extends from the San Luis Dam near Las Banos to a 
_ point near Kettleman City, Calif., where it joins state facilities to the south — 
_ (see Fig. 1). It is one of the major waterways in the western United States, 
more than 200 ft (61 m) wide, with water depths sup to 3: 33 fi ft ¢ (10 m a 


United States Geological in with the California 
' Department of Water Resources, began a study of subsiding areas in | California 
in 1956. Some of the findings concerning | shallow or near-surface subsidence _ 


_ Climate. .—One of the controlling factors | for the formation of collapsing soils 
is the climate. The southwestern San Joaquin Valley receives less than 7 in. 
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 *FIG. 2.—Typical Collapsing ‘Soil Mechanisms Which Hold Loose, Bulky odes i 


high-intensity storms, separated by long, dry periods (5), 
_ Geomorphology.—Shallow subsidence is confined mainly to the smaller coa- 
iL lescing alluvial fans that apron the foothills along the west | and south : side of | 
. the San Joaquin ‘Valley. Most of the drainage basins contain ‘unconsolidated te 
clayey sediments, have steep slopes, sparse vegetation, and short periods of — 
high rainfall. The intense rainstorms cause impact erosion, acquire a heavy 
load of loose, fine-grained material, and then move the material down confining, 


Pos steep-sided channels and deposit sediments 0 on the alluvial fans. This type of 


“deposits with a an open structure of high p porosity and low 

eu _ Texture.—Detailed studies of the texture of the collapsing soils revealed voids” 

_ between grains held in place by clay bonds, bubble cavities formed by air 
entrapment, interlaminar openings in thinly- laminated sediments, and voids left — 
_ by disintegration of entrapped vegetation. . Fig. 2 shows some typical mechanisms 

which hold grains together i in low-density soils. The predominant clay 


in all the subsiding fans is montmorillonite, which has a 
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— alinement of tl the San Luis Canal. These areas were oueay in the 
interfan regions between outwash streams and were considered to have been 
caused by the rapidly-deposited, low-density mudflows of fine- grained soils. 
"Since the San Luis Canal was to be constructed through these areas, it was 


is that the areas of shallow subsidence be delineated and the seriousness 
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FIG. 3.—Problem Areas of Land Subsidence and Route of San Luis Canal Through 


of the problem be evaluated so that pretreatment of the canal route could be 
= prior to construction. Approx 20 miles (32 km) of the canal passed 


through soils which required extensive investigation to determine the potential | 
gz shallow subsidence (see Fig. 3). Because of the size and importance of 
the canal, limitations in settlement were of essential consideration for the 
Tr maintenance of grade and the prevention of cracking in the linings and the 


e 


“the near-surface subsidence problem w was quite se serious. Several tests 
had been performed and some showed more than 10 ft (3 m) of near-surface — 
subsidence as a result of wetting. Much of the route of the canal was to pass 


type irr irrigation from 


| - Reconaissance surveys of land subsidence in the San Joaquin Valley showed __ 
| 
| 
4 SouTwERN \ PACIFIC eo; : 

| _ through areas where at least moderate, sprinkler- 


investigation procedures were used which wouid identif y the 
of those soils. These field investigations were supplemented with laboratory _ 
tests on representative samples taken at frequent depth intervals. Finally, several — 

field test ponds were installed to correlate the probable shallow subsidence — 
with the varying of ‘Seriousness by the explo ation program ani and 
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‘~Specific gravity = =2.70\ Lines 100 % 
Specific gravity =2. 60 


the liquid limit © 
for Evaluating L Looseness ant Prob noe 
undisturbed soil-sampling program was conducted with drill 
F- 100 ft (30 m) deep, spaced at 1,000-ft (305-m) intervals along the canal © 


-alinement through the subsidence : susceptible areas. Samples ' were taken at depth 
intervals of 5 ft (1.5 m) for the entire depth of the hole. Holes spaced at ae 
1-mile (1.6-km) intervals wre extended to the groundwater table, which was 
_ generally near 200 ft (61 m) in depth. Undisturbed samples were obtained at 
_ 5-ft (1.5-m) intervals for the entire depth of these holes also. For the San Luis 
Canal investigation, thin- d drive samplers v were used and provided 
undisturbed sample r recovery. ‘The method used was similar to that: described 


in a the American for Testing and 
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be us used to advance the drill holes because “water rand drilling fluid have 


Since different types of soil may be stable at quite different ranges of density, : 
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TEST DATA FOR ORLLL HOLE E-502 
EXAMPLE OF DATA FROM A DRILL HOLE SHOWING 
STABLE AND NON- SUBSIDING SOILS 
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DATA FOR DRILL HOLE E-514 


(b) EXAMPLE OF DATA FROM A DRILL HOLE SHOWING 
AND NEAR CRITICAL SOILS 


ORY DENSITY~pef SATURATION % LIMIT = 


Clayey ‘sone = 


-Fine, silty sand (sm) 
~Few grave! 
Silty sone (Su) 
Fine sand, cloyey (SC) 
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—Typical Examples ¢ of Exploration Data 
it was important to determine soil type bn to study the characteristics of plasticity, 
water-holding capacity, and the e effects of wetting in in regard to ‘overall canal 
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_ To explain and evaluate the critical condition for it was determined 


| 


a 
| 
~ 
| 
il 
of 
* 


and susceptible to sufficiently ¢ dense and unlikely to subside. 
One of the criteria for analysis, making use of the natural in- n-place density a 
and the standard liquid limit test, is slhown in Fig. 4 (3). Soils of different : 


= Case Il), saturation can easily cause a “liquid limit ahaeses, at which 
Zz the soil offers little resistance to deformation. If voids are less than this 
amount, as shown by Case III, the soil at saturation will remain in the plastic 
€ state and have greater resistance against particle shifting and only settle as 
a normal result of loading. If the voids are greater than this a amount, as shown 
by: Case I, saturation would ‘result in a moisture content in excess of the liquid © 
limit and the potential for collapse would be high. If collapse did not comm, 
- the soil would surely be in a very sensitive condition, - 
Using the aforementioned criterion of liquid limit and natural density, each 
_ drill hole was studied to determine if the natural densities were above or below 7 d 
i this limiting density, based on the liquid limit moisture content for each sample © 
tested. Although the soils must first be in a critically low- -density condition 
for a collapse of structure to occur, the existing degree of saturation was also” 
considered an po aU part of the analysis. Soils which were already at a 


Fis. 5 shows example graphs of the detailed data studied tor each drill 
_ First, the natural density was compared to the limiting density based on - 
liquid limit. Then, the natural moisture content was examined to evaluate the 
: existing degree of saturation. ‘Finally, the results of these tests were anal 
the descriptive log of the drill hole. 
Lasoratory Test Data 
| The logs of the drill holes were quite important for establishing confidence 
- in the laboratory test results, since there are practical limitations in the theoretical ‘ 
criteria. . First, it is assumed that the soils are uncemented, and second, | that 
they are sufficiently fine-grained so that the liquid limit test is applicable. To 
_ apply these criteria to cohesive soils containing gravel or coarse sands, considera: 
tion must be given to the fine soil fraction, because that is the portion to 
_ which the liquid limit test is applicable. Nonplastic sands are also questionable 
as to applicability of the criterion, unless” they" are fine — so that the 


s that soil density was of hi Stype of subsidence =f 
- ‘to occur, the soils must be sufficiently loose so that they are capable of — : 
when their particle-to-particle bond is weakened by wetting. Soils of different =f 
= require different quantities of water to reduce them to a weak plastic — 
¥ condition. This suggests that a simple criterion can be used to indicate the 
; limiting density at which any given soil will collapse upon wetting. Then, by 
! asucity OF water-noiding Capacity will Collapse a liferent densities. liqui 
| limit is a moisture content, determined by standard laboratory tests, which _ j 
| ~~ the weakest plastic condition of the soil or the condition at which ‘. q 
_ 6=s« itt is approaching the liquid state. When the soil has a low density such that 7 
q etting, whereas soils with low or moderate degrees of saturation (less than 
a 
| 


is not applicable | for soils with liquid limits less | than 20%. 
There are also practical limitations regarding the accuracy of natural density 
determinations. It is assumed that a small sample represents the average © overall — 
~ surrounding mass density. If frequent changes in soil layering or Sipe 
are present, then the natural density from relatively small samples may not 
fT s | 
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be representative. Also in such soils, small samples selected for liquid ae 
and specific gravity tests may not be truly representative of the entire sample. 
_ Therefore, numerous tests at frequent depth intervals were made for the San 
‘ Luis Canal investigations to establish confidence in average values. These factors 


— some of the principal causes of irregularities shown in some of the data, 


| 
ade 
bat 
| 
foo 
ved iso 
| 
at 
Or Of Consolidation lest Snaecimens Attar We Was bin lied 
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establish trends for the effect « loading o on 1 subsidence, laboratory 


in each drill hole, generally at 10 ft, 25 ft, 50 ft, and 100 ft(G@m,76m, © 
_ 15.2 m, and 30.5 m). The test loads were standardized so that samples from 24 
the 10-ft-25-ft (3-m-7.6-m), depths were tested under a loading of 10 Ib/sq 
in. (69 kPa) to represent an approximate overburden pressure. For the same 
reason, the deeper samples to 50-ft (15.2-m) depth were tested under 25 Ib/sq 
_ in. (172-kPa) pressures, and samples: from 50-ft-100-ft (15.2-m-30.5-m) depth — 
_ were tested under pressures of 50 Ib/sq in. (345 kPa). These test specimens __ 
| Were trimmed from the 3-in. diam (76-mm), thin- wall undisturbed yam tubes. 


_ belargely th that which would occur when any test specimen i is reloaded, priclay 7 


After final measurements were made under natural moisture niaiiies: the 


F settlement taking place under this condition was considered to ae collapse due 
_ to wetting. Many of the test specimens did not change in height, some showed 
‘further settlement due to wetting, and a few showed very slight expansion. “ale 
To compare the general results of these laboratory consolidation tests with — 
the aforementioned liquid limit-natural density criterion, the settlement observed 
when water was added is indicated on the chart shown in Fig. 6. Charts such 
as this were prepared for samples from shallow depths, intermediate depths, 
and from the 100-ft depth tested under the loadings described previously. ef 
__ These laboratory consolidation tests showed that the criterion of ages 4 density aa 
based on the liquid limit moisture content is 
soils where overburden pressures are small. At greater depths, mes as 100 
ft (30.5 m), these tests indicate that the limiting density should be considered 
- to be about 10 Ib/cu ft (160 kg/m’) higher than that indicated by the liquid A 
limit. - Generally, it it was found that soils most susceptible to were 


e less, ‘the criterion of limiting density based « on ‘the liquid limit is a ven very ‘useful F 
_ guide for assessing the =: of ¢ the soil and its probability of — 


= reaches of the canal to 5 be treated would be kept to a minimum but would — 

SY, include those areas which were critical. Where minor looseness was found, 
only moderate treatment might be required and areas which were me feed, 

fs stable would require no treatment. Using this approach, an attempt was made — 

: to determine for the San Luis Canal alinement areas that: (1) Were most critical; 3 

q (2) had critical soils extending to considerable depth; (3) had only near-surface 
critical soils; and (4) were previously delineated as having subsiding soils but | 
were investigation to be not critical b but have normal Properties. 
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all the results in a combined form. Fig. 7 shows the “method used to ohana 

the data from an overall examination of results for an entire suspected —. 
area. The charts are arranged to contrast the two most important features: _ 

(1) Whether in-place density is above or below limiting density, based on i 
liquid limit; and (2) 1 the degree of saturation | of the the soil at in- “place v water content. % 7 


density and the in-place degree of saturation is below 100%. Therefore, the 
«CANAL ALINEMENT--REACH 3 SUBSIDENCE AREA 


ESD Toto! irrigation about 0-9 aF/ac TEST POND 3. 
= Toto! irrigation over 20 aF/ac 


BETWEEN NATURAL DENSITY AND THE LIMITING DENSITY--P 


DEGREE OF SATURATION) 
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“Settlement meosur cane! route after 18 months of poraing 
«FIG. 7 —Settiement sila Ponding Compared to 
greater th the amount of black shading, the greater the likelihood of — 
and conversely, the greater the amount of gray shading, the better the . 
conditions. The estimated amounts of past irrigation, as well as the locations 
of drill holes along the alinement are also shown in Fig. 7. Only one small 
portion of the alinement is shown on the figure for illustrative purposes. The 
a group of 10 graphs on the left indicates that the soils from these holes would — B 
have very little subsidence potential. There is a prominence of gray area ng 
vy ry little black | shading showing. Densities app appear to be well above the | limiting 4 
— based on the nese limit. In some cases, Jess degree of saturation is 
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. concern. Soil conditions appear much more critical in the area ea shown on ‘the 

; is right in Fig. 7. Black areas are much more prominent and only a few holes 
show high degrees of saturation, 

_, It was observed that the most serious subsidence areas were generally | located 7 

ause 

_ where the estimated amount of prev previous irrigation was was small. Because irrigation 


Second embankment is to be constructed 
after ponding is completed 
Woter surfoce 


“SECTION A-A 
8. and Section of Test Ponds 


evaluating the seriousness of 
As mentioned previously, : several field test | ponds were constructed in areas. 
where had shown conditions most critical for ‘subsidence due to 


Bite 
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a 9.—Ponding Test in Progress Showing Settling of Banks 


_ wetting. wea — these ponding tests were used to ith the e exploration _ 
Program and laboratory tests. 
diagram of a typical test pond is shown in Fig. 8. As can be seen _ 
a the figure, each test pond had an embankment constructed on one side to represent — 
a typical canal embankment and to demonstrate the effects the embankment 
might have on subsidence. - After the ponding tests were completed, another . 
embankment was | ; constructed across the test pond, as shown by the dashed _ 
lines, to determine the amount of additional settlement due to the weight of 
the embankment. Each test pond had four gravel-filled infiltration wells to assist = 
in wetting at lower depths and to penetrate any impermeable soil layers. Each : 


test pond had a series of surface settlement platforms installed, both inside 
‘and outside the pool area (see Fig. 8). At least one deep settlement marker 


was installed at each pond to monitor subsidence at depth. Fig. 9 shows 
: panting test in progress and the associated subsidence due to wetting. 


Data from the test ponds had shown ‘ponding to be an effective method of 


octiageing and densifying these low- -density soils. It was therefore determined =a 


10. —Ponding Operetions to Subsoils Prior to o Construction of Canal 
yg ponding would be used prior to construction | to densify the soils along 
the canal alinement, through the areas most susceptible to collapse.  apbaider: 
- Ponding (as shown in Fig. 10) performed along the canal alinement aoa 
to construction caused subsidence of varying magnitude. _ Measurements or 
settlement devices spaced at. 500-ft (152-m) intervals permitted a direct comparison 
with the investigation data. The aforementioned critical area subsided 4 ft and + 
3 ft (i. 2 m and I. 5 m) in | many places, as shown a at the bottom of Fig. » would 


Critical areas. Stabilization by the ponding treatment of the loose subsoils down 
to the existing wi water ‘table near a of 160 1 ft (49 was 
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of water had been epplied dusing | perious irrigation 
operations. Collapsing of these loose soils was very effective in reducing the 
amount | of differential foundation movement after the embankments were 
constructed and the canal was operating. A total of about 20 miles (32 km) — 
i _ of canal alinement was treated by ponding prior to construction. __ — 
‘The mor movement of the phreatic su surface fr from the ponding was observed using - 
i a commercially- -available nuclear moisture meter for both the experimental plots _ 
oe and for those constructed during the actual ponding operations prior to construc- —_ 
tion (4). Casings in which the nuclear moisture meter was lowered to take 
_ measurements were installed to a depth of 50 ft (15.2 m) near test pond 3, . 
at distances of about 30 ft, 50 ft, and 80 ft (9.1 m, 15. 2 m, and 24.4 m) f from z a 
the pond; readings were taken for 17 months (see Fig. 11). Holes were > then 
drilled and undisturbed samples obtained to confirm that the "subsoils were 
saturated down to the water table at a depth of about 160 ft (48.8 m). satiate am 
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Range of | phreatic lines at preconstruction 

7 11.—Movement of Phreatic Line with Time Traced with Nuclear Moisture Probe : 

rivers 

the right-of-way line on the 100 

ft (30.5 m) from the ponds. Readings were taken after 3 months, 10 months, _ 


To take full advantage of the preconstruction 1 ponding operations, water is ; 


ponds. In the areas where the phreatic ‘surface 1 rose ‘pasting: 
terminated in as little as 1 yr, because the water table was approaching _ 
_ canal-bottom grade and no significant settlement was being observed. pan CR 
: _ A comparison of the phreatic lines after 3 months and 17 months of ponding — 


at test plot 3 and several locations the preconstruction ponds 


_ The San Luis Canal embraced more extensive and greater Hs Fle es of problems 


| 
preconstruction ponds were fille wi water, Casings tO receive 
_ the nuclear probe for moisture measurements were installed on 1-mile (1.6-km) | q ie 
> 
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- 40 ft (12.2 m) for fills and 50 ft (15.2 m) for cuts were required. In reach 
& expansive soils were encountered, and in reach 3, collapsing soils. It was 
anticipated that conventional slopes of 1-1 /2:1, that were used elsewhere, were 
steep for | these problem soils. Therefore, flatter slopes of of 2 or 2- 1/2: 


and this was also considered in the slope stability analyses. The canal was 
eventually constructed using 2:1 side slopes. 
For construction in the collapsible soil area where the subsoils were subjected = 
q preponding, special consideration was given to the probability of a soft, — 
saturated condition of the canal foundation. After ponding was terminated, a 
 dealanae period of several months was permitted to allow the initial collapse 
to be as complete as possible and to extend the effect to soil at greater depth. 
Also, these wet soils had to dry somewhat to permit workable conditions for — yr 


= slopes (4). The San Luis ‘Canal is “located in a a seismically- active area 
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“FIG. 12. genoa Results ene Test Section Using 50-Ton Pneumatie- Tired Roller 


& of initial collapse. Therefore, to - provide a firmer subbase to support the 
concrete lining, specifications called for rolling the bottom of the canal and» 
= the embankment foundations with eight passes of a 50-ton (445-kN), pneumatic- a 
tired roller. At a test area, the effects of this consolidation were investigated 
by sampling before and after rolling to determine changes in density. It was" = 
found that some increase in density occurred as deep as 3 ft (1 m), but that 7 
most densification took place in the upper 1-1/2 ft (0.5 m) (see Fig. 12). Rolling» 
of the canal side slopes was also considered desirable, but no feasible method 7 
was developed for accomplishing it. Since experimental data showed that 
embankments would continue to settle after they were built, _ specifications a 
a six-month period for stabilization after construction was complete 
_ before concrete lining was placed on them. In areas not considered sufficiently _ 
critical to warrant preponding, sprinkling was required to wet the foundation 
‘soils enough for settlement development in preparation for embankment con- = 


7 greater water depth, as much as 33 ft (10 m), and required considerable freeboard i . 
a 
a 


a The through which ‘the canal wa was ‘constructed have s soils 
. Thin, impervious layers are not ra 
uncommon. These prevent eatin and sprinkling methods (which collapse the aa 
soils) from being effective unless unfiltration wells are installed. The northern 
_§ miles (6.4 km) of the subsidence area were ponded without wells for more & 
than 1 yr. One bridge site has since - settled more than 3 ft (1 m), because = 
of an extensive impervious clay layer at a depth of 25 ft (7.6 m) which prevented © 
subsoil saturation (2). This is considered a minor problem in relation to the 
magnitude of the overall conditions encountered and the performance of the 
canal has been excellent. Because of deep subsidence due 


_ o structural damage ‘of the canal has been observed which could be attributed 
_ to further hydrocompaction of the near surface soils. The San Luis Canal was 
= during 1963-1968 and, as a result of the extensive soil investigations ms 
nd preconstruction treatment of foundation a soils, he has performed quite well. aE 


‘The: alinement of the San Luis Sonal in the San Joaquin Valley in California, “a 
a pie through several areas of low-density soils that are subject to collapse — 
upon wetting. These critical soils were delineated before construction by extensive i 
field exploration and laboratory | and field tests. Criteria f for ‘defining the limiting — 
- soil density at which collapse could occur were developed using the liquid d limit — 
_ moisture content and the in- place” density of the soil. These criteria were 
. _ Supplemented by a detailed study of drill hole logs, laboratory consolidation 
test data, and observation and measurement of settlement from test ponds. __ 
Nearly 20 miles (32 km) of collapsing soils along the alinement of the canal _ 
were ‘stabilized by ponding before | the canal was constructed. The surface soils ; 
"were generally saturated down to the water table which was at a depth ca 
about 160 ft (48.8 m). Further densification of the foundation soils was accom- ‘| 
= plished during construction by requiring eight passes of a 50-ton (445- bane 


pneumatic- tired roller on the canal invert and on the foundations for canal ~ 


- Preponding the ¢ collapsible soils was an effective method of stabilization. The 
‘ . canal was built during 1963-1968 and has performed quite e well to 10 date (1979). 
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ComPUTER SIMULATION OF CREEP OF CLay 


te ‘By Roland Pusch and Paul F 
dependent. 7 TAD :be == 
Creep of clay under undrained conditions may he of considerable importance 
_ practice, as illustrated by time-dependent settlement of foundations = 
_ overconsolidated clay, and “‘depth creep”’ in slopes excavated in a low- permeability, — 
4 Applied soil mechanics still makes use of empirical laws for iiss prediction _ 
of creep strain since no generally accepted creep theory is yet available. — - 
Object of this paper is to present an improved, physical creep model and its _ 
: mathematical analogy, and to examine the creep rate according to this model _ 
means of The yield 
to have a wide scope and it thus provides a promising basis for “representing 
“exep data, and a rational starting point for creep studies which involve structural 


‘Natural, soft, illitic pome consist of rather dense and rigid aggregates [Fig. 
ir 1(a)] which support each other, or are connected by links or bridges of small ' 
particles (2). Most of these are associated edge-to-edge or edge-to-face with Z 
direct mineral / mineral contacts, or separated by one or a few water molecule 
layers. The variable character of the c contacts ts implies : a spectrum of bond strengths, 
7 shear : stress is s applied, the induced internal stresses ‘will ‘shear ‘the weaker, 
interlinking particle groups, transforming them to ‘‘domains’’ with the particles _ 
= — 
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_ The interparticle distance in the domains tends towards « an eq 
in the stress field as determined by the local mechanical stress, the interparticle 
and molecular attraction, and the repulsion resulting from 
at small ‘interparticle. distances, or from electrical double- -layersa at larger distances. . 


of clay-adsorbed water, but also ‘ ‘aged”” domains remain the weakest 


structural members of the clay. "Displacements in the particle network in the | 

b course of creep may therefore be favored at these locations. Such comparatively — a 

j weak microvolumes may be considered as sustaining substantial deformation a 
and thus as main slip units in the evolution of Set their defo mations and 


AGGREGATE, 


FIG. 1.—Formation of Domains: (a) Natural, Undisturbed Condition: ade Local O 


displacements, at least partly, resulting from water molecule interactions in — 

form of breakage and reformation of hydregen bonds (3). 
a _ The weakness of the domains facilitates a relatively rapid initial deformation 
at the onset of creep which, however, becomes progressively retarded by several * 
| a mechanisms. Firstly, movements of the domains will be such that strong aggregates 
= and big grains will approach one another, make contact and strengthen the 
structure through mechanical interlocking. Secondly, strength may be gained a 


the reformation of water and the formation of new 


| 


of adjacent aggregates or grains has ceased, either because of interlocking : 
or pew they are temporarily not subjected to strain-inducing local stresses. 
- Thirdly, the displacements may lead to the formation of new domains through 
the yielding of bridging links which become ‘‘overstressed.’’ If the stress level 
is sufficiently high, such structural damage will accumulate and bulk a 
Wiiaid otni joval-» edt moit wolbuo wi sao 
The heterogeneity in stress and structure on the microscale will jointly result 7 
_ a distribution of heights of the energy barriers. The structural changes imply - x 
Z that the energy spectrum will will vary in the course of the creep. Thermal activation zz 


is assumed, for the creep p of soils. is known to be significantly temperature 
dependent. This suggests that, as in other materials where the creep is known > 
0 result from thermally activated glide, it is governed by a stress-dependent 
and time- “dependent ‘Spectrum of energy barriers, which control the | ‘unit slip” 


~ model of the creep of soils (4). In practice, the equations s obtained may not 
- be amenable to easy analytical treatment, particularly if one wishes to examine 
the influence of the energy-barrier distribution at the onset of the creep on 
the time dependence of the creep rate. A simple computer model has therefore 
2 devised (1), and its application to the creep of clay is here examined 4 P 
in the specific case when the initial distribution of energ energy barriers is “flat,” 


A basic assumption is that the dwell time, 6, of a slip unit at a barrier of - 
effective u, is given by the Arrhenius rate equation — 


— 


in u= the barrier height; = an atomic vibrational frequency of the 

order of 10'*/sec; k = Boltzmann’s constant; and T = the temperature, in 

At any temperature only a limited energy spectrum 

1 

will be of for with u < uy ‘the dwell- times will too short, 
‘for u > u, too long, to be ‘ ‘observable”’ ’ during the periods over which creep , 
tests generally extend in practice (4). 

ga If slip has been activated at a certain point in the clay, i.e., a barrier has _ 


been overcome, a contribution to the overall shear is made by | the i associated 


q 
| 1 | 


JANUA 


= be encountered red by the ‘same spreading slip-zone will be either higher or 
ete by an average amount 5u. Concerning the form of the spectrum, i.e., 
i the number n(u,t)3u of potential slip units per unit volume of material bed 
j up at time ¢ at barriers of height u, the spectrum is considered to be subdivided 
into equal intervals 8u. Now referring to Fig. 2, which shows consecutive a 
7 activation energy intervals, the time rate of change of the number n(u, Bu ¥ 
of flow units in the u-interval, i.e. , (an/at)du, is ‘seen to be determined on 
the one hand, by the outflow from the w-level into the higher energy-interval — i. 
a (1/2) n(u,t)8u v, exp (—u/kT) and, it is assumed, of the same number into 
the lower one, indicated respectively by the right and left lower arrows. my iil 


‘Tf all the separate contributions are added, one obtains 


nut u t) ex t) ex 

+ m(u — Bu, 1) exp | — 


al Assuming that 8u << u, — u,, the right- -hand side of the equation can be 


seen to tend to the secon pen differential of t) exp ( 


respect to u, as 8u — 0. 
= yields, on writing for n 


5 
If each activated is taken to make the same, mean, contribution 


stull Piast ad} = bee 


My ‘represents the upper limit of the operative sgeerem. In the integration it 
is assumed that the scheme of jumps shown in Fig. 2 is preserved at both 
ends of the cascade; thus u, is regarded as a ‘“‘generating barrier,’ the formation __ 
of new slip units being favored in regions where the activation energies for Ss, 
a ‘slip are low. These t new “units are taken as ones the otherwise absent net 


= 
uy, will occur with ‘such a low frequency that their contribution to the creep : 


: "over the period of observation may | be neglected. One may regard uy as an 


kil) with 
ar 
4 
| 
| 
| 
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ae 


“‘absorbing barrier.”’ These considerations lead to the result (4) that the wie 
probable value of u in the spectrum will be equal to about 25 kT, i.e., about 


A solution of Eq. 3 which le leads to creep laws geneeelly observed in solid aa 


vplt 
= of +t.) 


in orw the appropriate constant of proportionality depends on the deviator 


_ The evolution of creep according to this creep model is here examined by a 
% means of a computer code (1) for numerical il integration of | the equations governing 
the thermal activation, using discrete energy intervals in the thermally activated . 
movement of independent slip units across energy barriers. 
_ Code.—The code is based on the assumption that the probability that a slip — 
unit body have a kinetic energy sufficiently high to overcome a barrier encountered, 
e~“/*” per attempt. If the attempt frequency of the slip units is then 
number of encounters in ‘time At will be v pAt. 
_ The probability, p, that a given “unit held up at a barrier at time t, is still © 
neld up by the barrier time + Atistherefore 


—u/kTs 


If a large number of independent slip ‘units are contained in independent — 


potential wells of depth u, , the which ‘has crossed the 


_ We assume that the range of barrier heights u, — u, increases from uy ef 


steps to u>. . With a total of N steps we therefore have 5u = 
= u,)/(N oo 1). Activated jumps will” be | assumed to take a given s slip unit, 
>. held up at a barrier of height u, with equal probability (1/2) to a barrier of 
height u + 8u or u — Su. This single-step process characterizes the numerical-in- > 
4 tegration scheme which yields the distribution law, n(u,t), for the barrier heights. 
= In the present work we shall assume that n(u,0) = const for all u-values - 


j the range. This boundary condition will underestimate the n-population for low 


3 u- 1-values compared with the distribution given by E Eq. 5, as is s is readily ‘inferred = 
on considering n(u, 0) i in both cases. od} sroods or 


| 
| 
— 
43 
| 
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The integration scheme for the model to the cascade 
used in the derivation of 3i in the foregoing, is as follows 
4 


_ assuming in the present instance that there is no net influx from outside the am 

range, i.e., from Bas. This implies that no units are generated at weak 


(low u). Furthermore _ bie sels 
2(t- + At) =p3(t) - frp 
Py(t + At) alt) (14) 


condition will be. unimportant in practice, for will be chosen to 
_ imply dwell times long compared with the period over which the creep is being 


Parameter 


considered. 


Time step for integration, in 0<DT<o 


DT: 
-TNDS: Total number of slip units 0 < <0 
TPRINT: Time step for printing results _ <TPRINT. 


u, = Minimum energy, in electron volts — 

= Maximum energy, in electron volts — 
7 TEMP: Temperature, 7, in degrees Celsius EM 
BC: Boltzmanns’s constant, eV/K 


‘Attempt frequency, < FREQ < 


As we were interested principally in the shape of the creep curve it 
t 


possible to choose the values of all pecnenetns appearing only in wale constan 


7 
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FIG. 3. of Benter per du tntorwel of Energy 
: Spectrum i in Course of Creep (0.1 eV-0.9 eV) : 


FIG. 4.—Total Number of Slip Units Held up at Spectrum Intervals at Different Times — 
after of Creep: 10° sec; 5 x 10° sec; (3) 10° sec; (4) 5 x 10* sec; 


— 

| 


of (e.g. see Eq. 7) arbitrarily, ‘for conv convenience in computing. 
‘The values used were as follows: NEB = 10, DT = 10, TMAX = 3-10° = 7 


= = 100 sec, FREQ = 1/sec, TEMP = 293 K, TNDS = 1,000. i” 


first attempt was to investigate a a rather spectrum: 0.1 eV 
.: = 0.9 ev, ‘i.e., a range well covering that of the hydrogen bond. Fig. 

4 shows the successive redistribution of slip units held up at the lowest boundary 
and at u,, u,, u, and u,, which represent the boundaries between 


lower four energy intervals of the spectrum. Since 10 intervals were taken, 
there were 100 slip \ units held up at per ati = $e 


FIG. 5.—Total Number of Slip Units Held = a in Each Spectrum Interval 
at Different Times after Onset of Creep (1) 10° sec; (2) ‘10° sec; (3) 5 x 10° sec; 
(4) 10° sec; (5) 2 x 10° sec (0.1 eV- 


- successively. It is obvious, however, that intervals higher than 0.6 eV-0.7 eV 
_ do not participate significantly with this u-spectrum at room temperature, which ey 
: indicates that only a narrow part of the spectrum has to be considered . This 

was confirmed by other runs using 0. l eV-l. 5 eV and 0.1 eV- -0.6 eV spans, ria 


the latter being shown by Fig. 5. qu biel mes 


_ As to the creep rate, it is assumed in the theory that every jump makes 


the same contribution to the strain. The creep rate is thus obtained from the 
values | of nu, (Figs. 4 and 5) using | the relation uae 


| 
all 
| 


for various values of 
Fig. 6, showing the —_ based on Fig. 5, indicates that the computed relation 
is of the type given by Eq. 7. 
The relatively high ¢ ,-value, i.e., the comparatively low | initial creep rate, =) 
and results from the assumed “‘flat’ * shape of the n(u,o) spectrum, points to 


- the importance of the initial u-distribution in determining the form of the aaa : 


_ The simple model used in the present computation ‘yields creep” curves 4 
“ee, also ‘reflect the “typical behavior of m many illitic clays 2 as observed in — 
— 


* 

vot Wedron = 


4 FIG. 6. —Relationship between ya = ten and time, for 0.1 eV-0.6 eV 
(Straight tne is Obtained on with t, = 7,000 sec instead of ¢) 


_ The microstructure of fine-grained soils is highly heterogeneous, implying ri 
large spatial variation in the strength of the more or less continuous particle — 
network. In creep induced by a constant deviatoric stress, the number and 
_ heights of thermally activable barriers to slip will thus vary in space and time. 
The resulting evolutionary change in the distribution of the spectrum of barrier 
# heights is encompassed by a stochastic model of the creep of clay previously | 
_ described by the writers. A simple version of this model was used in a computer — 


ee to determine the form of the creep curve at constant stress and = 


— COMPUTER SIMULATION — 
| 
q 
; 
isin many cases rather small (5). Further use of the code with different 
rs 
Jogarithmic, in good agreement with experiments. The use of various boundary, 


conditions, particularly defining the barrier spectrum at the onset of creep, | i 
should facilitate the diagnostic use of the computer program in studies: of the =. 


Miles, G., and Felthamn, P., ““-DENEB—A Computer Code for Numerical Integration 

4 of the Equations Governing the Thermal Activation of Discrete, Independent, Dislocation © 
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The following are need this paper: 


shear iis in contribution due to single ac soit jump; 


fraction of slip u units eee crossed barriers of height u; 


Boltzmann’s constant; 
ay ; 


energy it interval step; 
number of flow units or barriers; 
probability; 
abbreviation; | 


y 


= dwell time of slip u unit; 

v = slip frequency; and ion 

= atomic vibrational f requency. dor 
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‘STRESSES | IN SOIL AROUND Yarra 7 
COMPRESSIBLE PILES 
S. Sankaran,’ N.R. Krishnaswamy,” and B. | Chandra’ 
An attempt has been made to evaluate the stresses the 
shear distribution using linear elastic theory as suggested by Mattes and Poulos — 
(3). The influence of the compressibility and L/d ratio of the pile on the stress" 
_ distribution have been studied for »p = 0. 5, to show the difference inten ¢> 
a very flexible and a rigid pile. The side shear distribution thus obtained for _ 
a particular ratio of L/d and compressibility of the pile is divided into uniform | 


and linear variations (), and is finally summed up, to get the net stresses. = 
This is a step-by- -step procedure, and the double integration of the Mindlin’s: ’ 
- equation. (4) involved is solved partly analytically and partly numerically in q 
_ acomputer, using Gauss Quadrature formulas. The analysis is confined to vertica on 
stress, although other stresses could be evaluated by the same technique. ti oe es: 
The analysis is more realistic than the earlier methods (dd, 2) f for the following | 


Feasons: (1) There is no assumption: regarding the side shear distribution; 
the length to diameter ratio of the pile is taken into account; (3) the compressibility 
7 of the pile, se peernan by the stiffness factor, A (ratio of Young’s modulus 
of the pile, E,, to the Young’s modulus of soil, E,), for a solid pile, is taken 
_ into account; and (4) further, the stress at the tip of the pile > can be be evaluated. 
The assumptions made in the analysis are: (1) The pile base is assumed 
a 4 be a smooth disc, across which the base load is uniformly distributed; and 7 
o 2) the disturbance of the continuity of the elastic halfspace due to the presenc a 
‘The table of stress coefficients is presented to cover a range of stiffness 
“factors representing a highly flexible to a very rigid pile. In tables Lag a 
stresses can be interpolated for any intermediate flexibility of the pile. The 
tables are presented for selected ratios of m (ratio of depth to the length of 
the pile, z/L) and n (ratio of radial distance to the length of the pile, r/L). s 
The L/d ratio and stiffness factor are varied from 10-100, and 50-50,000, 
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0’; and 5 x 10*, respectively. 


= =8=©6TABLE 1.—Vertical Stress Coefficients (x 1,000) for L/d = 10 (+ Tension, Otherwise ‘etal 
| 
or (2373 1,855 | 1,435 | 863 | 519 | 91] 49] 17] 7 | 
$23 | 474 | 365 | 270 |153 | 97] 66] 28) 12) 
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244 | 237] 219 | 197] 175] 124] 85 | 39] 18 
| 249 | 230 | 207} 185] 130) 89 | 41 
268 | 266 | 264 | 251} 232 | 209) 187} 132) 89) 42) 19 
127 127 | 126 | 125 | 122} 116] 109] 100) 77| 57] 30] 1s 
216 188 | 181] 171 | 159] 145} 111] 81 | 41 | 20 
199 198 197 | 194 130 180 | 167] 153 M16 | 
29} | 200] 199 | 196 | 192] 182| 169] 154] 117] 85 | 43 
Note: The four values noted one below the other correspond to \ = 5 x 10'; 5 x 10°; 5 x | 


TECHNICAL NOTES 


at 


n, Otherwise 


Compression) 


m| 00 018 | 02 | 03 | 15 
) (6) | (7 (8) (13) (14) 
5, | 217 | 64 
126 
123 
386 


242 
192 


+1,583 
+1,656 
+1,€71 
+647 
108 
615 
(1,626 
37,769 | 6,375 
45,045 | 7,741 
7,917 
3,392 
4,104 
4,196 
1,685 
2,044 
2,090 
210 


ar 


Note: The four values noted one below the other correspond to A = 5 x 10'; 5 x 107; 5 x ul 


“4 Coctficionts (x 1.000) for 25 (+ 
d 
7 
540 480] 395 332 179] 133] 63] 28) 1 
| ‘18 173 | 143 | 74 & 7) 
om +132 | +34] 11s | 159 | 186 |172] 144] 76 6) 7 
0.6 332: | 401 | 387 | 306 145] 4] 
+303: | +137 +17) «145 | 235 | 284] 254] 204] 102] 49} 11) 
+657 | +407 | +227, | | 268 [262] 221] 118] 13] 
+704 | +442 | +254) | 152 | 266} 263| 223] 120] 59} 13) 
| +479 | +187 96 198 | 224} 189] 146] 68] 32) 8| 2 
+599 | +78| 355 | 455 | 407/318] 244/125] 65/17} 5 
+568 | +6 | 456 | 353 271 | 142] 75 | 20) 5 
+569 +1 | 461 |357| 274] 144] 76] 22] 5 
09 414 70 199 | 194] 167] 135] 71] 35} 10) 
937 | 1,028] 839 | 658 | 445 247/131] 70/20) 6 
1,373 | 1,375] 1,054 | 799 | 516 |370|278|145| so} 24| 7 
1,421 | 1,415 | 525 |375| 282] 150] 82] 24] 7 
1.0 «447 201 | 168} 145] 122] 70] 4 
3,328 | 2,095| 1,133 | 762 | 460 |327| 247| 134) 75] 23) 7 
4,117 | 2,617| 1,416 | 940 | 548 283] 152] 85/27) 9 
4215 | 2,681| 1,450 | 961 | 559 |386| 287| 154] 87} 27] 9 
| 555] | 164] 133] 112] 68] 
2,735 | 2,147| 1,337 | 896 | 503 |341| 252] 137] 79/26] 9 
3,320 | 2,616 1,634 1,097 | 607 | 402| 293| 157} 90/30] 10 
lt 3,395 | 2,675| 1,675 | 1,121 | 619 | 410| 298] 159] 91 | 31 | 10 
379) 300 238 | 164 | 128] 105] 66) 40 | 14 | 
a... 1,559 | 1,410] 1,107 | 852 | 526 | 357] 261| 29] 10° 
1,893 | 1,713 | 1,346 | 1,036 | 635 | 425] 307| 163] 95] 34] 12 
a 1,935 | 1,752] 1,376 | 1,059 | 648 | 433] 313] 165] 96/34) 12 ff 
205 | 200! 186 | 170 | 139/114) 94] 62) 40) 16) 7 
gag | 735 718 | 696] 640 | 574 | 445 | 339] 261/147) 89/34) 16 
993 | 890 870 | 774 | 694 | 536 | 406| 311 | 171] 103 | 40] 16 
1,014 909 888 861} | 547 | 415] 317] 175} 104 | 40 
| | | 423] 387 | 335 | 282] 145] 92] 39 17° 
$m | 527 | sar] 466 | 403 | 338| 279] 171] 108] 45] 19 
«$82: | 538} 533 | 525] 503 | 476 | 411 | 344] 284] 174] 110] 46 | 20 
on 9S 95 | 94) 92 | 9] 84] 77| 69] SI] 37} 19] 9 
| 298 | 296 | 294] 287 | 277 | 253 | 226] 197] 136] 92] 42] 20 
«382: | 358 | «356 | 345 | 333 | 304 | 270] 236] 161] 108 | 49] 21 
390 366. | 364 | 361] 352 | 340 | 310| 275| 164] 110] so| 
72) | 67] 62] 58] 46) 35) 19] 10 
| 932 220 | 219 | 210 | 197] 181] 164] 122] 88] 4a] 22 
| 264 | 263 | 261 236 | 217] 196| 145] 109] 25 
269 | 268 | 266] 262 | 256 | 240] 221] 200] 148] 106] 26 


TABLE 3. —Vertical Stress Coetticients (x 1 000) for L/d = 50 (+ Tension, Otherw se 


+339 
(+441 
+1,524 +713 
+1,897 | +818 
+835 
265 | 
1,021 
1,148 
2,134 
3,980 
4,278 


0 
AS 


Seaat 


0.2 7,025 | 5,055 | 3,627 | 1,812 | 
2,987 | 2,201 | 1,649 | 933} $30 | 182 | 73] 36] 12] S| 
O42 544 376 | 256 | 129 | 78) 53) 25) - 
430) 422 | 386 | 294] 214 120 | 56) 27) 12) 
1,400 | 1,302 1,041 | 779 | 391 | 185] 87] 17] 
1,564] 1,268 | 1,055 | 756 | 554 | 309 | 183} 114] 43] 19] 4 | fh 
| «201 | «216 | «218 | 201 171 | 138} 71) 33] 
#80] 56 126} 163 | 183 | 169) 142] 76] 36) 
24 
. 
+198 | 364] 504 | 460 | 361/279) 147) 77/21 | 6 
naa +31 | 62} 99 | 123 | 116] 99] 27) 7) 20 
506 | 444 | 335 | 258 | 104] 55 | 15 
0 1,201 | 987] 765 | SOS | 365 | 276| 146] 79 | 23 
x 1,315 | 1,062 | 816 | 531 | 381 | 287] 153] 83 | 25 
8,542 134 | 89| 9 | 98] 96] 85] 52] 28 
56,007 1,301 | 704} 493 | 326 | 246/193] 106| 58 | 18 
93,948 2,480 | 1,331 | 890 | 529 | 371/278] 150| 84) 27 | 8 1 
99,440 2,670 | 1,432 | 954 | S61 | 390/291] 157) 88] 28 | 9 
507 | 348) 275 | 212 | 137] 105 88 | 82) 74) 49) 299) 10) 4 
(3,510 | 2,447) 1,940 1,496 | 906 | 602 | 349 | 247 | 190] 107 | 61 | 20 
6,237 | 4,313 | 3,443 | 2,675 | 1,632 | 1,075 | 589 | 392 | 287/155] 89/30 | 10 
6,506 | 4,610) 3,683 | 2,863 | 1,748 | 1,150 | 627 | 415 | 302] 162] 93| 31 | 
24 172 128] 107 8S | 74) 67) 46) 29) | 4 
4,215 1,211] 1,114 592 | 367 | 255] 192] 108] 64] 22 
2,564 | 2,138 | 1,969 | 1,771 | 1,372 | 1,041 | 627 | 417/301] 160] 93/33 | 120 
(2,652 | 2,286 | 2,105 | 1,894 | 1,467 | 1,113 | 668 | 443/318] 168| 98] 35/13 
(101 9 | 92] 87) 82 | 72) 63) 
961 522] S10 | 494 | 453 407 315 | 242 | 189] 110] 67 | 26 
1,016 911} 889 | 860 | 787} 703 | $37 | 404} 308} 169] 101) 39 | 16 
1,042 919 840 750 | $73 | 430 | 327} 179 | 107 | 41 | 16 
A | | | 52] 47] 36] 26) 13) 6 
308 | 303 | 291) 275 | 239 | 202| 169] 107| 69 29 13 
| 534) 528 520 | 498 i 404 | 357| 278] 169/106} 44) 19 
| 569} 563 554] S31} SOL | 431 | 359) 295] 179] 112} 47 | 200 
«50 49 | 49] 48 4 | 43 32} 24); 13 7 
4 213} 212 | 210 | 205] 198 | 182 | 162/143] 99] 68} 32] 15 =f 
| 361} 359 | 356 | 347 305 | 270 | 235| 160] 107 | 48 | 22 
(399 385| 382 | 379 | 369] 356 | 324 | 287| 169] 113 | 
| 0) 3] 39 | 35] 33) 27) 22) 13 | 
“2 263 is8} 157 | 156 | 154] 150 | 142] 131} 119] 89] 65/33} 17 J 
“a 265| 264 252 236 | 217 | 196] 144] 103 | 51 | 25_ 
(291 280 | 279 | 274] 268 | 251 | 230] 208] 154| 109| 54] 26 
Thee four values noted one below the other correspond to = 5 As 
10°; and 5 x respectively. 


TABLE 4. ._—Vertical Stress Coefficients (x 1,000) for L/d= = = 100 (+ Tension, Otherwise 


4,988 3,693 | 1,860 
3,946 | 2,850 | 1,489 


191,043 
194,408 


Rone ws 


4 


4 Note: The four values 
10°, and5x 10‘,respectively, 


i 
5,539 773 | 210 | 58| 16} 10 
1,484] 1,161 | 916 569 355 | 150 | 77 | 47 19 ‘Si 
«453: 331 234 | 125] 77) S54] 25] | 
2,657} 1,849 | 1,916 | 1,859 900 | 386 | 194} 22] 7 
2,232 | 1,857 | 1,562 | 1,111 785 | 382]185| 21] 7 = 
aa] | 538 | 428 351 | 2451177] 130] 61 | 27) 
97 | 133 | 196 | 195 | 171} 140] 73 | 34 
06] | 228] 313 (389 | 391 379 | 263 | 128] 45] 17 | 
194] 303 312 | 299 | 209) 139) 49] 19 
#202] +73 | | 159 | 237 | 282 | 252 | 203) 102 | 49 
+541) +337 +186 | | 160 | 264 | 260 | 220) 118] 58 
| 16S | 166 | 142] 113] 55] 25 
+1,089 +578 | +269] 48 | 169 | 211 | 181] 141] 65] 30 
1,675 | +752 | +208 | 294 431 | 406 | 321 | 247| 126 | 65 
| +1,812] +790 | +188 499 | 456 | 358 | 276] 144] 75 | 20] 
09 103 | 131 | 136 | 121] 101] ss] 28] 8| 
| +905| 228 300 37 | 181: | 158] 129] 67] 34] 
| 4653} 693 931 | 817 655 | 448 | 330| 250) 132] 70/20] 
| 1,249 | 1,032 | 799 | 524 | 377 | 284] 150] 81 | 24] 
1.0 318 186 | 145] 125 | 120 | 114] 103] 89| sa] 30] 
46,894 1,547 379 | 211 | 173 | 153 | 136] 115] 67] 36] 
3,246 | 2,011 | 1,090 | 743 | 458 | 328 | 249] 135] 75 | 23 
4,076 | 2,541 | 1,373 | 922 | 549 | 384) 287) 155] 27) 
11] 256 | 186 | 160] 128 | 113 | 100] 90] 79] 30] 10] 
725 | $47 | 324 | 221 | 149] 124] 105] 65] 37) 12] 
4,984 2,915 | 2,252 | 1,365 900 | 502 | 341 | 252] 138] 79] 26] 
6,118 | 4,524] 3,600 | 2,788 1,694 1,114 | 610 | 406 | 297) 159] 91 | 31] 
| 138] 133 | 126] | 104 | 90) 80) 71) 49) 30) 
of 1,025 460 294 228 | 153 | 119 98} 63] 38/13] S| 
2,046 1,815] 1,670 | 1,500 1,159 | 879 | 533 | 359 | 263] 142] 83 | 29 
2,512 (2,288 2,060 | 1,850 | 1,430 | 1,083 | 650 | 432 | 312] 165] 96 | 34) 
14] 88 | 86 Ss 82 79] 72] 65) 59] 43) 29] 
i 206] 202] 196 | 181 165 | 133] 108] 89) 58) 38] 16] 6 
813 755 | 731 | 669 | 598 | 458 | 266] 148] 39) 35] 14 
998 | 927 | 897 820 | 559 | 420] 319] 175 | 104] 40| 16 
1264} 121 103} 90] 78] $4) 37) 17) 8 
46 346 | 289 | 239] 147| 93 | 39 
574 556] 550 | 541 518 | 489 | 420] 288] 175] 110] 46] 200 
wat ‘et al 46] 43] 40) 33] 25] 14] 7 
294 | 90] 86} 80] 73| 66] 49] 35] 18 
316 | 309] 307 | 304 296 286 | 261 | 232} 202] 138] 93] 42) 200 
373 | 360 | 348 | 317| 280] 244] 165 | 110 | So 3 
#4 4. 40 39 | 39 38] 36] 34) 29) 23] 13 
lan | wl 6 | 67 | sol ss] a3] as) 9 
al 26 204 216 | 202 | 186 | 168] 125} 99 | 44 2 
272 | 268 | 261 | 245 | 225 | 203] 150] 106 | 52 6 


Lenten,’ England, Vol. 16, No. 3, 1966, pp. 231-255. 
< 2. Geddes, J. D., ‘“‘Boussinesq Based Approximations to the Vertical Stresses Caused 
by Pile Type Subsurface ice Loadings,”” London, England, Vol. 19, No. 
3. Mattes, N. S., and Poulos, H. G., “‘Settlement of Single Compressible Pile,”’ Journal _ 
of the Soil Mechanics and Foundations Division, ASCE, Vol. 95, No. SMI, Proc. 
-Mindlin, R. D., ‘Force at a Point in | the a “4 "Physics, 
Vol. 7, May, 193 1936, pp. 195 202, 
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r Discussions may be submitted on any Proceedings paper or technical note published in any 
Journal or on any paper presented at any Specialty Conference or other meeting, the Proceedings __ 
of which have been published by ASCE. Discussion of a paper/technical note is open to 
_anyone who has significant comments or questions regarding the content of the paper/technical — 
oe Discussions are accepted for a period of 4 months following the date of publication . 
of a paper/ technical note and they should be sent to the Manager of Technical and Professional — 
‘Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period 1 my 
be extended by a written request from a discusser. 
‘The original and three copies of the Discussion should be submitted on 8-1/2-in. (220-mm) | 
by Il-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of 
a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 

_ of manuscript including figures and tables); the editors will delete matter extraneous to the — 

- subject under discussion. If a Discussion is over two pages long it will be returned for shortening. a 

7 All Discussions will be reviewed by the editors and the Division’s or Council’s Publications _ 

_ Committees. In some cases, Discussions will be returned to discussers for rewriting, or they” a 
may be encouraged to submit a paper or technical note rather than a Discussion ve es i. (4, a 

Pi Standards for Discussions are the same as those for Proceedings Papers. A Discussion a 
subject to rejection if it contains matter readily found elsewhere, advocates special interests, 

r carelessly prepared, controverts established fact, is purely speculative, introduces personalities, 2 
or is foreign to the purposes of the Society. All Discussions should be written in the third ee 
_ person, and the discusser should use the term “‘the writer’’ when referring to himself. The : 
author of the original paper/technical note is referredto as “the author.” 

_ Discussions have a specific format. The title of the original paper/technical note appears — 

at the top of the first page with a superscript that corresponds to a footnote indicating the 

month, year, author(s), and number of the original paper/technical note. The discusser’s full a “a 
— should be indicated below the title (see Discussions herein as an example) together with Pan 

his ASCE membership grade (if applicable). 

The discusser’s title, company affiliation, and business address should appear on the first 
g of the manuscript, along with the ngerenny paper number of the original paper / technical = 


Note that the discusser’s identification footnote should follow consecutively from the original 
note. If the paper / technical note under discussion contained footnote numbers 


¥ - appeared in the original paper/technical note use the same number used in the original. al 
It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 


the Publications of ASC ASCE for more detailed information on preparation and bon 


ge ne Figures supplied by the discusser should be designated by letters, starting with A. This also = 


REME- VaLue Prosiems or Limrrinc 


Discussion by A. Barron,’ F - ASCE 

_ This paper presents an analysis method for stability problems at limiting 
hee we using the method of armneranen calculus. This calculus i is not new, 


had they been included, \ would have improved the ‘clarity and 


1. No list of assumptions is given. Near the end of the introduction the 
authors state that none are used. Certainly the use of the effective stress concept _ 
as expressed by Eq. | is one. Another is the use of safety factor, F, to apply 
not | only” for the overall stability, but also as a point safety factor when F 


_ 2. No illustrative examples are given. To aid the comprehension by most 
engineers a detailed example would be most useful. 

_ 3. No comparison with results of others are given. Near the end of the section 
“Slope Stability Problems” the authors state the Tesults of their method 


__ 4, No discussion is included indicating the influence of the soil stress- strain 
properties on the validity of the method. In fact, in the second paragraph of 
the “Summary : and Conclusions” the authors State that the method * ‘takes into 

and pore-water 


# Nowhere in the paper are the strength parameters defined. The writer assumes 
since an effective stress equation is used (Eq. 1) that the strength parameters 
F are those obtained from consolidated-drained (S) shear tests. The writer concurs | 
4 in the application of Eq. 1 to failure conditions where F is equal to unity; 

but for cases where F is greater than unity the drained shear strength for the 
- nonfailure effective normal stress is generally not available. See Barron (10,11) 
2 Gould (14) and Johnson (12,13). Except for dry, cohesionless soils the potential 

shear strength of a wet soil not at a failure condition is the undrained —_— 4 


tests are, therefore, | ‘the consolidated-undrained (R) or “un- 


- drained (Q) as the conditions to be analyzed indicate. In the writer’s opinion, 


+ | TOeiaber 1979, by Michael Garber and Rafael Baker (Proc. Paper 14901). eye 
* Consulting Engr., 62 Horshoe Rd., Guilford, Conn. 06437. 


y the profession, additional shou given. it 
| 
| 7 
q 


is the most misused vee soil mechanics. oP 
test is made, the effective stress path f for a | dilative soil has a positive slope — 
and at failure the pore pressure can be negative (gage) ‘resulting in a shear 
strength at failure larger than that given by Eq. | with the effective ioe 
7 stress that which exists on the potential failure surface at the initiative of the — oa 
est. Because large ‘negative pore pressures (gage) can occur in a test using © oe 
back pressure to prevent pore-water cavitation, the full shear 
given by this test may not be conservative tot use ina stability study. 


If the soil is compressive (strain-softening), then the effective stress . 
a rises to the failure strength curve with a negative slope. Large positive pore © 


pressures can exist and the shear strength at failure is less than that indicated - 


‘ed! 
In addition, if F, the : average safety factor, is larger than unity, then the # 
“stress conditions in the soil mass are highly indeterminate and the concept kee 
of a limiting equilibrium condition with a reduced shear strength being applied a 
to a given problem is certainly an assumption, but one which is commonly _ 4 
_ made without, however, claims being made for any greater degree of ecco, 
 Anexcellent discussion of the difficulties of analyzing slopes is given by Johnson _ 
(12,13). It would be desirable if the authors define Dirac’s function and Lagrange’s <- 
In conclusion, the writer questions the advisability of publishing theoretical a 
_ papers of this kind without it being so written that it can | be a 
a large portion of the s membership. 
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vig by Ryszard J. 


The authors applied a a ‘limit equilibrium approach to the | problems 
of slope stability and problems of limiting magnitude, and location or direction 


of a load and external moment. The most important result is a basic theorem 


“Adjunct, “Inst. of Geot technics, 
5 


of Wroclaw, ul. Glogowska 9/72, 


of iit which states that the critical value of the 
; parameter (F = the | factor of safety; P = an external load; M = an = 
% ‘moment, etc.) is independent of the normal stress distribution along the critical g 
Slip line. Simultaneously, the variational analysis establishes the existence of _ 2 
_ two (rotational and translational) modes of failure mechanism. _In the case of . 
_ the rotational mode of failure the potential slip lines are a log-spiral but if 
the translational mode of f failure exists, the potential slip. lines are straight lines. a 
analysis” was performed without any a priori assumptions regarding slip 
line and normal stress distribution; all the equilibrium equations are satisfied aa 
. no aditional assumptions were in the investigation process. The results “ 
_ stated in the paper are very interesting but are not surprising, = = 
_ So, approximate methods for solving the problems of soil stability may | be 
divided int into (15, 16, 18): (1) Methods of limit analysis based on two limit theorems; 
and (2) methods of approximate satisfying of conditions of equilibrium and — 
yield criterion only in the definite points or regions of the material considered. | 
These method are called the limit equilibrium methods. 
In the limit analysis for elastic or rigid-plastic material obeying the associated 
flow law the validity of two basic limit theorems is proved, making use of | 
the principle of maximum “plastic work. According - to the statement of this 
oa in order to properly bound the “‘true’ " solution, it is necessary to 
_ find a kinematically admissible failure mechanism (velocity or flow field) in 


order to obtain an upper bound solution. A statically admissible stress field, 
Satisfying all the equilibrium conditions and nowhere violating the yield criterion, 

will be required for a lower bound solution. If the upper and lower bounds - 

coincide, the exact value of ‘the collapse, or limit, load (factor of safety, etc. * e 


An upper bound solution may be obtained (see. Ref. 18) by: ray Comparing 
_ the sum of the total work of external forces and total work of body forces 
with the total internal dissipation of energy (this method is called herein energy 
approach); or (2) making use of equilibrium conditions of the field of forces 


(stresses) associated with the assumed kinematically admissible failure mechanism .. 


_ Both the approaches are equivalent. It can be observed, that the equations 
of the work balance, may be treated as an equation of virtual work and, thus, 
aS an equation that expresses the condition of global equilibrium. A number 
of examples solved by means” of the two approaches are discussed in | Ref. 

The upper bound ‘equilibrium approach is more advantageous in ‘relation to 
_ the analysis starting from work balance. Namely, it provides more information 
about the stress distribution inside the material or the force (load) distribution 
along contact lines. There are, however, difficulties in applying this approach 

in t the case layered soil profile and _ composed mochenioms, | which 
consist of rigid zones (blocks) and soft zones. 

_ In the case of associated flow law the energy ps within the failure 
@ a mechanism is independent of the stress distribution. Therefore the determined _ 
value of limit load (and all other parameters) will be independent of the internal 
stress distribution. In the kinematically admissible failure mechanism, the potential — 
slip lines are either log-spirals (rotational mode) or nes (translational 


¥ 

| 


The aforementioned properties of the upper bound solutions and kinematically 


_ admissible mechanisms follow from the associated flow law. The same properties 
_ are obtained by the authors by applying, in the light of the assumed classification — ; 
of methods, the upper bound generalized equilibrium approach. Their considera- tee 
tions, however, are limited to the case of failure mechanisms consisting of 
rigid block which undergoes either rotation or translation, = 
| The upper bound energy approach allows the consideration of complex failure 
_ mechanisms composed of a number of rigid blocks which undergo either 
translation | or rotation, as well as of mechanisms containing soft zones, which _ 
are deformed during plastic strain. Some examples on composed failure mecha- 

nisms have been given by Gudehus (17), Chen (15), the writer and Mroz (18), 

3 To conclude we may add that conventional limit equilibrium methods og 
on static equilibrium of forces (or stresses) acting on the yield region are in 
- many cases equivalent to the upper bound method (see equilibrium approach). 

_ However, the assumed failure mechanism must be kinematically admissible and — 

_ it must be associated by the appropriate volumetric dilatation. When the failure : 
mechanism is not kinematically admisible from viewpoint of limit analysis, 7 
_ may be neither upper bound not the lower one. However, each upper bound a 
_ (equilibrium approach) solution is always equivalent to the limit it equilibrium ‘j 
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Discussion A A. Lucefio® and E. . Castillo® 


The authors state a value pri of LE as a search 
- for a pair of functions y(x) and o(x) that realize the extremum X,, of X subject 
to the satisfaction of the three LE equations (3). sate va 


_ However it can be easily demonstrated that functional X [y(x), o(x)] is, in 4 
general, Therefore there exists a potential sliding line y(x) and 


a normal str tribution | a(x) giving a a safety factor as small a as. desired, A 


*Escvela Técnica Superior de Ingenieros de Caminos, Canales, y Puertos, opto. de 
mame Aplicadas a la Ingenieria, Universidad de Santander, Avda. de los Castros, 

*Escvela Técnica Superice de Ingenieros de Caminos, Canales, y Puertos, opto. de 


— Aplicadas a la Ingenieria, Avda. de los Castros, s/n—Santander, Espafia. - 
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E _ In order to give an example of the unbounded character of functional K 


ag X [y(x),o(x)] let us consider Eqs. 3a, 3b and 3c for the case of a ne in 


7 


in which the following pase parameters have ‘been used: N = c/y H, 
S = o/y H, = tg(o), X = x/H, Y = y/H, Y = y/H, Y’ = dY/dX, 

_ where X,, X, = the abscissas of the end points + the = sliding line; and ae 
= the height | of the slope. 

By making X = (Baker and Garber, 1978), 29 g gives: 


according to Garber problem can be stated as minimizing 
32 subject to conditions described by Eqs. 30and31. 
_ One way of demonstrating that Eq. 32, subject to Eqs. 30° and i31,d does 
not attain an absolute minimum is by selecting a sliding line Y*(X), not ‘y 
_log-spiral, and by using the Ritz method to check that the functional F* [S(X)] 
by substituting the equation of this sliding line i in Eq. 32 and Eqs. 
30 and 31 has no absolute minimum. 
_ As it is well known, Ritz’s method consists of a discretization of the ae 
: ae of possible solutions. So, a system of basic functions [/,(X)] is — 


According to this assumption the [S(X)] becomes a function 
of n-2 variables from which its minimum must be investigated. ree 
effect, let us consider a slope by and 
geotechnical characteristics: c = | Tn/m* 3 ) = ‘Tn/m’; 


Ak The characteristics of the slope, in adimensional form are: N = ati 


0, 


— ‘ 
| 
| 


the following potenti slip line: = 3x? 1/3 and the system 


of basic functions: /,(X) = x’; l (X) = =X; 1,(X) = 1 are selected, functional 4 . 


F* [(S(X)}], according to Eqs. 32, 30 and 31, becomes a function of one ne variable 


Ja, +@- 12F*)a, + (18 18F*)a5 = 


Thi his system can be solved in a,, Se 


to —1/3. This shows" that Eq. 32 subject to | in Eqs. 30 


+ + + +0, 080 + 1,080 


Nevertheless, it is worthwhile to make some comments: 
1. Though the safety factor must be positive, this fact is not cndliaiite aitced 5 
‘ie the statement of the problem, and as a consequence the system, Eqs. ae 
(36 and can n be ‘solved for negative values of F*. 


is not as a ‘consequence the s: system of Eqs. 35, 3 36, ond 37 
solutions not satisfying this condition. 


_ One way of taking into account these contains in te ent of 
alt 
= 


Eq. 41 | defines functional K? = F which must be minimized subject to the 


conditions described by Eqs. 42 and 43. 
Mia extremals of this new Problem are are exactly the same as those 4 the 


(37) 
d 
| 
a 


 fanctional does not attain a relative minimum either (the first writer, 1979; 
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dine 


iscussers alps interest in the p paper. 
The discussion by A. Luceno and E. Castillo suggests that in the extreme- value cl 


- exist a pair of fanctions y(x) and o(x) * “giving (a) safety factor as small ys 
desired, and in consequence the extreme problem is incorrectly stated. is To 
ee it, the discussers provide a a‘ “counter r example” that, in their understanding, 
clearly disqualifies the writers’ presentation. example refers to a slope 
according to the writers’ solution, has a safety factor of 1.8. The discussers 
Bese at random a potential slip line y(x) and assumed the pressure distribution _ 
to be presented by a polynomial with three unknown parameters 
as. Substituting their functions y(x) and o(x) into > the writers’ equilibrium 
: equations the discussers derived a system of three linear equations with respect © 
to three unknowns a,, a,, and a,. This system can be solved for almost any 
value of the safety factor, F. It seems therefore, that the discussers can always 
point out a pair, y(x) and that not only gives Fless than 1.8 but, moreover, 
The point is, however, that there are natural restrictions on the value of 
F as well as on the character of the a(x) function, and the discussers re: recognize ‘ 
_ this as they make an attempt to modify the problem presentation. These restrictions — 
are: (1) F = O (the safety factor is non-negative by definition); and (2) biz) 


while o # (a condition supplementing the Mohr-Coulomb 


‘discussers. know, ‘their incorporation does not change ‘the variational ‘solution. 
If we consider the discussers ‘‘counter example”’ in view of these if ) 
a ’Sr. Research Assoc., Transport Tech. Lab., Dept. of Mech. and Aeronautical Engrg., § 
OR Room 203 C. J. Mackenzie Building, Carleton Univ., Ottawa, Canada KISSB6. a 


Sr. Research Fellow, Soil Mechanics Dept., Faculty of Civ. Engrg., Technion-Israel 


= 


q 
a 
7 
Hiem OF LF nresented Pv the writers the extremization parameter ttunctionai 


e 
pressure distribution S(x) assumed by the -discussers was in 
spected in the range 0 F < 4 [Figs. 5(b) and 5(c)] against restriction 
/ ag _ which in nondimensional terms is rewritten as S(x) = N/w. It was found that ad 
_ the minimal S always occurs at the end points x, and x, [Figs. 5(a) and 5(5)] 
Aa of the slip line y(x). As it can be seen from Figs. 5(b) and 5(c), S(x) Satisfies © 
ss restriction 2 for F = 4 but fails to satisfy it for any 0 = F < 4. This implies i 
+: F=4 is the minimum the discussers can claim for their example without 
violating: the ‘obvious. restrictions previously "mentioned. As 4 is greater 
1.8 the ‘‘counter example’ * clearly fails. Possibly in expectation of this the _ 
J discussers refer to another example that is supposed to demonstrate their point. — ties 


The is to two papers in Spanish one in » English 


0.4; 


FIG. §.—Investigation of Discusser's Exampl Discusser's Potential ‘Stip Line: be 


(Symposium in India). The writers’ efforts to these popes have been 
umsuccessful and thus comment is reserved. 
Va _ The variational solutions proved to satisfy restrictions | (which | in general a 
ee has to” be stated as X = 0, in which X = the extremization parameter) and im 
ian 2 as demonstrated by solutions presented by the writers (5). Taking the opportuni- 
om the writers would like to introduce an additional (kinematic) restriction that ha 
= be imposed on the shape of y(x). This third restriction states that the 
= is possible only when the curvature of y(x) decreases in the direction > 
of sliding. Thus » there ar are “proper” and “improper’” ‘potential s! slip lines, and 
“3 used by the discussers [see Fig. S(a)) obviously belongs to the class of 
“improper.” * The variational (extremal) Y(x) always satisfies “wi; third restriction. 
Indeed, using Eq. 85 the curvature p of Y(x)i is expressed as 


As: was 2s shown in the paper, the radius vector r increases in the direction a 


= (see Geometrical Analysis section), consequently p, (Eq. 44), decreases i 


The discussion by R. A. Barron questions th the manner in which the material “4 
: is presented. The writers readily agree with a number of points raised. The 
| So text was about 25% longer and it included more detail explanations i. 
_ about the essence of the method and its limitations. We had, however, to follow — 
reviewers request to considerably shorten the paper. ay 
The discussion by R. J. Izbicki provides a profound analysis 0 of the writers — 
method the of view of the limit The writers 


Page 1158, paragraph 1, line 3: Sh Should res read ‘ ‘the iigahe x direction,” ‘instead 

Page 1164, line 10: Should Tead “(rotational mode of failure),’” instead of 
mode of failure,’ 


DISCUSSION, 
ing corrections should be made to the original paper: 
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